Historic, Archive Document 


Do not assume content reflects current 
scientific knowledge, policies, or practices. 





Reserve 
aSB953 
. 1477 
1989 


Forest Service 
Equipment 
Development 
Center 


Missoula, MT 





February 1989 
3400—Forest Pest Management 
MTDC 89-2 


~ User Manual Extension 
For The Computer Code 
AGDISP MOD 5.0 










AGDISP 








United States 
Department of 
Agriculture 


NATIONAL 
AGRICULTURAL 
LIBRARY 


Advancing Access to 
Global Information for 
Agriculture 





eer rereeatipeneeienee cent vanie nesters TT UE nreumcmniermnmmemeiianiemeremmmmeemnmrmmemnrememremrssieseo 
> 3 


User Manual Extension 
For The Computer Code 
AGDISP MOD 5.0 






AGDISP 






i a tifa ap tp od 


— Nae 
om ae am «am am a= == —_ om =< 
“fe 





Prepared by 
Milton E. Teske 


Continuum Dynamics, Inc. 
P.O. Box 3073 
Princeton, NJ 08540 


Prepared for 

U.S. Department of Agriculture 
Forest Service 

Technology and Development Center 
Missoula, MT 59801 


Robert B. Ekblad 
Project Leader 
January 1989 


(This work was done under contract 53-0347-00952 by 
Continuum Dynamics, Inc.) 


Foreword 


This report is published as a part of the USDA Forest 
Service program to improve aerial application of pesticides, 
specifically by using pesticides and delivery systems 
tailored to the forest environment. The program is 
conducted jointly by the Technology and Development 
Center, Missoula, MT, and the Forest Pest Management 
Staff, Washington Office at Davis, CA, under the 
sponsorship of State and Private Forestry. The MOD 5.0 
version was cosponsored by the U.S. Army Dugway Proving 
Grounds and U.S. Army Chemical Research and 
Development Center, Research Directorate, Physics 
Division. 


Details of the aerial application improvement program are 
agence in two Fores Seis upeports, 





Planning and Analysis System (Forest Pest ‘Management 
Rpt. FPM 82-2, February 1982, Davis, CA). 


A system of computer models has been developed to 
optimize spray program design and operation and assess 
environmental risk posed by aerial spray operations. 
AGDISP is one of these models, and this user’s manual 
covers the most recent version of AGDISP, MOD 5.0. 


AGDISP was developed under the leadership of the National 
Aeronautics and Space Administration with support by the 
Forest Service to include special forest spray conditions. 


When NASA's Agricultural Aviation Program was terminated, 
the Forest Service continued to support and direct further 
development and improvements of the AGDISP model. 


The changes and additions that have been incorporated into 
MOD 5.0 are listed on page ii. 


The MOD 5.0 is operational on the USDA Forest Service 
Data General MV-15000 and is installed with GKS graphics 
and the IBM PC/XT/AT personal computers. 


SUMMARY 


The AGDISP computer code predicts the motion of agricultural material 
released from aircraft, including the mean position of the material and the 
position variance about the mean as a result of turbulent fluctuations. 
Developed under sponsorship by NASA, U.S. Department of Agriculture (Forest 
Service) and U.S. Army Dugway Proving Grounds, this code operates efficiently 
and is user-friendly, with many of its features validated against wind-tunnel 
and flight test data. This document is the Mod 5.0 User Manual to the AGDISP 
code, and includes a description of its current enhancements and operating 
instructions. This version of the AGDISP code is operational on the following 
computers: Data General (with GKS graphics); Digital Equipment Corporation VAX 
(with Calcomp and Tektronix 4025); and IBM PC/XT/AT (with proprietary plotting 
programs for dot matrix printers). 
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MODIFICATION LEVEL 


AGDISP code development has lead to the following versions of the code. 


Mod. 


Mod. 


Mod. 


2.0 


3.0 


4.0 


Operational on a Control Data CYBER 175 at NASA-Langley. 


Includes: all of the basic program development: fully 
rolled-up vortices or Betz roll-up, simplified 
terrain modeling, WAKE plot file entry, models 


for propeller, helicopter, crosswind, 
superequilibrium turbulence, canopy, vortex 
penetration into canopy, and material 


evaporation. Also includes a stand-alone program 
(AGLINE) to construct the equivalent Gaussian 
distribution. 


Graphics: Tektronix 4025 and 401X terminals using NASA- 
Langley graphics software calls. 


Operational on a Univac 1108 at Ft. Collins (U.S. Forest 
Service). 


Improvements: helicopter modeling with transition to rolled-up 
vortices; AGLINE calculations as a menu option in 
AGPLOT. 


Additions: discrete crosswind velocity profile; nonzero 
deposition height; composite deposition plots (to 
16 plot files); canopy penetration by helicopter 
downwash; plot option to plot material diameter 
time history. 


Graphics: DISSPLA at Fte Collins, with appropriate 
subroutine calls for Tektronix terminals. 


Operational on a VAX 11/785 at Dugway Proving Grounds (U.S. 
Army). 


Improvements: Betz roll-up procedure and _ propeller’ model 
revised; equivalent Gaussian distribution 
selection criterion is a_ program’ decision 
involving material vertical velocity. 


Additions: models for wide body effects, simple vortex 
circulation decay, jet engines, multiple 
powerplants and parameterized evaporation; 


default input file option; plot option to plot 
material vertical velocity time history. 


Graphics: CALCOMP at Dugway, with appropriate subroutine 
calls for the pen plotter. 


ii 


Mod 


5.0 


Operational on -a Data General MV-15000 at Missoula, and on IBM 
PC/XT/AT personal computers. 


Improvements: 


Additions: 


Graphics: 


Revised solution procedure eliminates integration 
stepsize dependence on material decay time 
constant. : 


Axial variation added to all models; ground 
Sprayer; continuous deposition; canopy 


deposition; deposition on objects. 


GKS at Missoula; dot matrix printer output on IBM 
PCLXTAAT 


pin let 
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NOMENCLATURE 


plant area fraction 

concentration 

effective canopy drag coefficient 
material or aircraft drag coefficient 
maximum body diameter 

material diameter 

target significant dimension 

inertial impaction efficiency 
gravity 

canopy height 

inertial impaction parameter, Eq. (48) 
material mass 

ground deposition 

mean square turbulence level, equal to <uu> + <vv> + <ww> 
radius 

Reynolds number, Eq. (8) 

aircraft semispan 

wing planform area 

time 

aircraft thrust 

fluctuating local fluid velocity 
mean local fluid velocity 

aircraft flight speed 

incremental axial velocity 

friction velocity 


fluctuating material velocity (u,v,w) 
x 


air 


Pair 


NOMENCLATURE (Cont'd) 


mean material velocity (U,V,W) 
relative velocity, equal to [Uy ss V; | 
swirl velocity 

downwash velocity 

aircraft weight 

fluctuating material position (x,y,z) - 
mean material position (X,Y,Z) 
canopy displacement thickness 
surface roughness 

reference height 

leaf capture efficiency 
circulation 

circulation at the wing centerline 
Kroneker delta function 

wet bulb temperature depression 
von Karman constant 

macroscale length 

forward advance ratio 

kinematic viscosity of air 

density of air 

density of material 

material standard deviation 
evaporation time scale 

material relaxation time scale 
turbulent time scale 


inertial impaction second parameter, Eq. (49) 
xi 


‘ 1. INTRODUCTION 


The original development of AGDISP was motivated by a desire to determine 
how aircraft-unique wake and propulsion characteristics affect the ground 
deposition pattern of aerially released material. Material motion was 
computed by a Lagrangian formulation of the trajectory paths and the solution 
of ordinary differential equations. Equations for the ensemble averaged 
turbulent variance of the released material were developed using a locally 
isotropic spectral density function. Idealized flow field models were 
constructed for aircraft vortex wakes, helicopter flow fields, propeller 
swirl, terrain, crosswind, plant canopy and superequilibrium turbulence. The 
result of this work was the Mod. 2.0 version of AGDISP as reported in Ref. 1 
and 2. 


Enhancements to this code include additional models for wide body effects, 
jet engines and vortex circulation decay. Refinements to the code, noteably 
to the helicopter model and the graphical presentation of the results, 
proceeded as the code was made operational at additional computer sites (Mod. 
pe enGtee 3? MOd 34.0. “Ret. 4). This document contains the latest model 
developments compiled in Mod 5.0, operational on the Data General computer at 
Missoula, the Digital Equipment VAX at Dugway and IBM PC/XT/AT. 


This manual is organized as follows. Section 2 summarizes the equation 
development and the models currently available in the AGDISP code. Section 3 
discusses the operation of AGDISP and includes a listing of the code 
restrictions, error and warning messages that may arise. Section 4 explains 
the input card options to AGDISP, and Section 5 explains the operation of the 
companion plotting package AGPLOT. Section 6 examines the test cases used to 
exercise the AGDISP models, while Section 7 reviews the usage of AGDISP and 
AGPLOT on the present computer systems. Section 8 highlights the selection of 
background turbulence level. The Appendices contain an explanation of each 
subroutine in the code, flow charts of code structure and a listing of 
significant variables. 
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2. SUMMARY OF MODELS 


A Lagrangian approach is used to develop the equations of motion of 
discrete material released from aircraft, with the resulting set of ordinary 
differential equations solved exactly from step to step. The formulation of 
the material equations, particularly the ensemble averaged approach used to 
develop turbulent properties of the material-atmospheric interaction, is 
detailed in Ref. 1. This section of the User Manual examines those variables 
and nomenclature that are a part of the AGDISP code results. 


The aircraft is assumed to be in level flight near the surface, releasing 
material into a (y,z) plane normal to the flight direction, -x , as shown in 
Figure 2-1. The AGDISP code solves for the transport of this material within 
this plane, until it deposits upon the surface or is carried aloft by wind or 
vortex motions. Accuracy of predictions requires accurate description of 
background environment (mean winds and turbulence). The need for accuracy 
must, however, be balanced by the requirement of computational speed and ease 
of calculation. Simulations need to be performed simply and repeatedly under 
a variety of background conditions, with minimum computer demands and setup 
requirements. When the need is justified, the solution for the background 
environment may be found by a larger, more detailed code, such as the NASA 
WAKE code discussed in Ref. 5. The mean velocity and turbulence profiles 
predicted by these codes may then be read by AGDISP to predict released 
material behavior. As an alternate approach the AGDISP code has been fitted 
with extensive simplified flow field options to permit the prediction of 
material motion in a wide range of idealized background conditions. These 
models are summarized in this section of the User Manual. 


The released material is assumed to be spherically shaped. The material 
flight path as a function of time after release is computed as the locations 
(X,Y,Z) for all material included in the simulation. The material velocity 
is denoted by CURVEW) F The material is affected by the background 
environment; however, the environment is assumed to be unaffected by the 
presence of the material. The interaction of the material with the turbulence 
in the environment creates turbulent correlation functions for the position 
and velocity, <yv> and <zw>.; for the velocity variance, <vv> and <Cww> ; 
and for the position variance, <yy> and <zz> . The square root of these 
last two variables gives the horizontal and vertical standard deviations of 
the material motion about the mean described by Y and Z. 


Material Equations of Motion 


In the AGDISP code the released material is assumed to be influenced by 
two forces: 


a) weight, and 
b) aerodynamic drag. 


The governing equations are derived in Ref. 1 and repeated here for 
completeness. 
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EES a C1) O=< (5) % Kae eG and Uj are the ensemble averaged ith 
component of material position, velocity and local fluid velocity, 
respectively, while Xj if v; and u, ‘are the fluctuating ith component of 


material position, velocity and local fluid velocity, respectively. 


Inherent in the material equations is a relaxation time T which is 
essentially the e-folding time for the released material to come up to speed 
with the local fluid velocity (for Vy to approach and equal U;)- The 
relaxation time is defined as 


Dp 
na BS Dee. (6) 


Eee 
3 Gp par wake 


P 


where D is the diameter of the material, Pp is its density, P.;, is the 
density of air, V,,) is the relative veloctty between the material and the 
local background, and C is the material drag coefficient. The semi- 
empirical formula for Cy (Ref. 6) is 


oo 


24 0.63 


—4.. 1.38 
Caleta (1 e201 I 7Re 4 22h UieRex: Gu) (7) 
D Re 
where Re is the Reynolds number defined as 
DV l 
Re oe ole ASSES (8) 
v 
air 
The material diameter Dp is affected by evaporation (Ref. 7) as 
dD oa 
ances ° 
p ; arate tery, 
e 
wherevpt. 8 the e-folding evaporation time 
eel ee 10 
“e ~ BAO (10) 
where AQ is the wet bulb depression in °C and 8 is defined as 
g = 84.76 [1.00% 0.27 Ret! >] (11) 


in units of m2/sec cat Guar 


In a recent set of tests (Ref. 8), the evaporation model described by Eq. 
(9) was compared with wind tunnel results. A typical comparison is shown in 
Figure 2-2, illustrating favorable agreement with the test data. 


Equations (1) - (5) cannot be solved without specifying relationships for 
the quantities <xju4? and <uj;vj> , the correlations of the material 
position and velocity fluctuations, respectively, with the local fluid 
velocity fluctuation. These expressions are developed by integrating their 
ensemble averaged frequency spectra using a spectral density function for 


transverse velocity fluctuations in isotropic turbulence (Ref. 9). The 
results (detailed in Ref. 1) give 
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where tT, gis the travel time of the material through a turbulent eddy of 
scale A , adjusted for the passive tracer limit 


+E — (15) 
wel - 8 4 


The quantity q2 


is the mean square turbulence level in the fluid. 

With the position and velocity information available for all released 
material at any time during the simulation, Eqs. (1) - (5) may be integrated 
exactly for the solution at the next time step. This step result may then 
serve as initial conditions for the next step integration, etc., through the 
entire simulation. The assumption that the background conditions Ue 
<x uy? and <u,V4? are constant across a time step permits solution of Eqs. 
ab - (5) unrestricted by the size of the released material. Computational 
solution times are reduced significantly, and heretofore unsolveable material 
diameters (including the passive tracer limit, Dp = 0 ) may be tracked by the 
equations. 


Flow Field Modeling 


The behavior of the released material is intimately connected to the local 
background mean velocity U, and turbulence field q2 through which the 
material is transported. The AGDISP code has been configured to accept user 
specified mean velocity and turbulence flow fields in neutral atmospheres, but 
also contains a number of simplified models for the flow field velocity and 
turbulence levels behind aircraft. These models are summarized below. 
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FIXED-WING FULLY ROLLED-UP TIP VORTICES 


When an aircraft flies at a constant altitude, the aerodynamic lift 
generated by the lifting surfaces of the aircraft equals the aircraft 
weight. The majority of the lift is carried by the wings, and generates one 
or more pairs of swirling masses of air (vortices) downstream of the 
aircraft. If the roll-up of this trailing vorticity can be approximated as 
occurring immediately downstream of the wing, then the mean velocity field 
that results may be simply characterized by the aircraft semispan, s 
circulation, [ , and load distribution (Figure 2-3), 


bd 


FIXED-WING BETZ ROLL-UP 


‘When the wing loading cannot be approximated by a single vortex pair, or 
the aircraft is flying sufficiently fast that roll-up cannot be assumed to 
occur immediately downstream of the wing, then the Betz methodology (Ref. 10) 
can be employed. The solution procedure (documented in Ref. 11) relates the 
swirling velocity distribution in the vortices to the details of the wing 
spanwise load distribution by assuming that angular momentum is approximately 
conserved. For complicated wing planform shapes, a vortex lattice analysis 
must be invoked to generate the wing spanwise load distribution (Ref. 12). In 
a simulation involving the Betz methodology, the vortex sheet rolls into the 
vortex as a function of time. Thus, the unrolled sheet will contribute to the 
ambient fluid velocity during the roll-up process. 


VORTEX CIRCULATION DECAY 


The flow field generated by a simply loaded wing includes four vortices: 
the two from either wing tip, and image vortices below the surface (Figure 2- 
4). The local ambient velocity is the vector sum of the four swirl velocity 
components, each written as 


V = pel. (16) 
Ss Tits 


with its direction measured perpendicular to a line connecting each vortex 
centroid and the observation point. The observation point could be a material 
position or the location of another vortex (which will also move because of 
the other vortices present in the flow). An inviscid vortex of constant 
positive circulation strength fT. = PF will descend toward the surface while 
increasing the separation distance between its companion vortex. 


In the atmosphere, turbulence will act to decay the vortex strength. A 


simple decay model, developed in Ref. 13, results in a vortex decay of the 
form 
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Figure 2-3. 


Swirl velocity distribution of a fully rolled-up wake 
as a function of spanwise load distribution and as 
measured from the center of the resulting vortex. 
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T = iy exp(- eats (17) 


For vortices out of ground effect, a typical value for A = 0-41 ; while a 
detailed examination of Program WIND data in ground effect (Ref. 14) gives 
Aq = 0.56 m/s . 


HELICOPTER IN FORWARD ADVANCE 


The helicopter model includes both the hover downwash and the tip vortex 
pair by partitioning the helicopter weight between the two effects as a 
function of time. The hover downwash model is taken from actuator disk theory 
for a propeller and may be written as 


ae, 
= 18 
FW. 279 348 Wy (18) 


while the strength of the vortex pair may be found from 


Gi Py Ns .= 20 e RUT (19) 


where W is the helicopter weight, Wq is the downwash at the rotor plane, 
R is the radius of the rotor and [I is the tip circulation strength. The 
function F is determined by the expression 


kur ot | 
F = exp {- 5 ) (20) 
oTR 


where wu is the forward advance ratio, oO is the solidity and I is the 
solution to Eq. (19) for F=0.. The constant k_ relates helicopter roll-up 
(around a circumference at the blade tip) to fixed-wing roll-up. Predicted 
results from a helicopter wake model (Ref. 15) suggest that the downwash flow 


field quickly transitions (within two blade revolutions) into a vortex pair. 
A typical value for k=7. 





The tip vortices influence the flow field everywhere, while the downwash 
affects only the spreading region beneath the rotor and the fluid column 
directly above the rotor. Within the boundaries of the rotor blade and the 
dividing streamline (Figure 2-5), the velocity is assumed to follow a 
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| Rotor Downwash | 
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Vortex D 


Vortex 
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Figure 2-5. Schematic of the helicopter flow field model. 


stagnation point flow, 


while the vortex pair created by the helicopter is 


assumed to follow the dividing streamline. 


Material released ahead of the helicopter (spray boom forward) is assumed 
to encounter a streamline pattern similar to flow around a circular cylinder 


(the actuator disk downwash, 





R? IR y" 
LRU Det ee 
E r 
2U_Roxy 
V=- ; 
Cc 


Ref. 16) 


(21) 


(22) 


until the- material crosses the plane of the helicopter shaft centerline. 


JET ENGINE 


The jet engine exhaust is modeled as a turbulent circular jet using the 


similarity analysis reported in Ref. 1/7. The results give the axial and 
radial velocities as 
3 Vv T/e os, 1 
Yaxial ~ On Ex 2 2 (23) 
Cl +en /4) 
Vv 14/3 wPair Sule n/t) (24) 
i 
radial 4 T x ie: Soyo es 
with 
Beles SN: 
4 Vin € x (25) 
6 = 0-0161-/T/p., (26) 
where T is the jet thrust and x is the downstream distance (corrected for 


the virtual origin). 


turbulence level 


The radial velocity 
fluid velocity, while the axial velocity 


Vradia] becomes part of the local 
Usxia] iS also used to evaluate the 


q wx 0.2034 aera Cos 


where the proportionality constant is determined from the centerline decay of 
turbulence in a free jet (Ref. 18). 


PROPELLER 


The propeller is modeled as an actuator disk, where the incremental 
velocity AU over the flight speed Us is related to the actual thrust 
produced by the propeller 


THs 2 WoT =RoAUCU,, + AU) (28) 


if 


where R is the radius of the propeller. In steady flight the thrust equals 
the aircraft drag, so that 


] 2 
T= Cy Gs Pe eas: (29) 
where S is the wing planform area and C is the aircraft drag 


coefficient. Combining Eqs. (28) and (29) to eliminate thrust, we obtain 





(30) 


An alternate expression for € , defined in Eq. (26), may be obtained from 
Ref. 17 as 


Combining this equation with Eq. (26) gives an effective thrust level for the 
propeller of 


yoy) 
= BZ 
Tere P5575 eds AU (32) 


to permit the use of Eqs. (23) - (27) for the propeller when Eq. (30) is 
substituted. The swirl velocity generated by the propeller is assumed to be 
linear in Rout to its maximum value, obtained by integrating the axial flux 
of angular momentum, and then zero for larger values of ‘r . The resulting 
integration yields 


VEL por Ea ee a4) 
PNR (U,, + Av) 


where n is the propeller efficiency and 2 is the propeller rotational 
speed. 


WIDE BODY EFFECT 


The presence of a significant fuselage shape can generate additional mean 
velocity and turbulence in the vicinity of the aircraft. For simplicity, we 
assume that the fuselage cross-sectional area can be modeled as_ an 
axisymmetric body shape B(x) , where x is measured from the nose to the 
tail. Slender body theory (Ref. 19) then gives the additional radial velocity 
contribution from the changing area of 


Ue ; 
= Bex) (34) 


Vv ‘ ~ 
radial 27Yr 
where the prime denotes a derivative in x . 
The turbulence level, consistent with the axisymmetric wake created by the 


passage of the body, is found by matching the data in Ref. 20 with a 
turbulence model (Ref. 21) and assuming an exponential decay 





2 

2 2 1 r 
q- =#)q exp|[- 5 ( ) (35) 

body 2 R/2 

where 
1.43 
2 2. 7X 

rc 0.003U" bas ) (36) 


and 


(37) 


is the wake radius to half-width, D is equal to the maximum diameter of the 
body and x is corrected for the virtual origin. 


MEAN CROSSWIND 


In a neutral atmospheric surface layer the horizontal velocity follows a 
logarithmic profile 


gn (z/z.) 


gn (2/20) eee 


V(z) = V(z.) 


where vCz) is the known horizontal velocity at a given altitude z, , and 
z is the surface roughness. Surface roughness is typically taken to be 
1730th of the actual physical roughness height. With a linear integral scale 
of turbulence (A = 0.65z) , the turbulence level becomes 


2 
V(z,) 


Cr 0.845 lin Ce_/z,)! 


(39) 


SUPEREQUILIBRIUM TURBULENCE 


The detailed effects of turbulence are obtained by invoking super- 
equilibrium turbulent transport theory (Ref. 22). Superequilibrium refers to 
the second-order closure turbulence transport model limit where the velocity 
correlations are able to track their equilibrium values. The nonlinear 
equations to be solved for the turbulence may be written in index form as 


dU 6. 
m 4 bie 929 fated) 69 ez 
aX [Spee aoe oA <uyu,? | z A [<uju,> 3 oe 1 A A LD 


sl j 


Since the local mean flow gradients are known, the system of equations 
represented by Eq. (40) may be solved exactly for <uju;? to determine q* . 


PLANT CANOPY 


The canopy flow field model is an approximation of the second-order 
closure turbulence model discussed in Ref. 23. For a canopy of height he 
and plant area fraction A(z) , the mean crosswind velocity above the canopy 
is 


gn Cz/ zi) 
Cc 


(41) 
Qn (z/z.) 


V(z) = V(z) 


wheréyszeojis the displacement thickness of the canopy. Within the canopy the 
mean velocity and root mean square turbulence levels are assumed linear with 
height, in good agreement with data (Figure 2-6). The presence of the canopy 
alters the downwash stagnation flow below the helicopter by moving the 
effective-surface up to Zz, . Vortices entering the canopy will have decaying 
circulation strength of the form 


fisge 1 
a) MG ye: 
1 + —— f(t) 

21S 


similar to the atmospheric decay effect in Eq. (17). Here the function f(t) 
integrates the cumulative effects of ACz) onthe .circulation of the vortex 


t h 
Co e 
to) ag J J A(z)f ,dzdt (43) 
6 H -Ah 
Cc 


where Ce is the canopy drag coefficient, Ah is the penetration depth of 
the vortex and f, is the fraction of the vortex within the canopy. 


TERRAIN 


Surface slope is modeled by assuming that the ground plane remains 
straight but may incline. Any vortices in the flow field will have a modified 
image vortex system to maintain zero flow through the ground. The helicopter 
downwash model and the position of the dividing streamline are also altered. 
Crosswind and canopy effects are assumed to remain parallel to the surface. 


GROUND SPRAYER 


A ground sprayer is modeled by eliminating the presence of the aircraft 
vortices and being very specific about the direction and initial velocity of 
each nozzle. 
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Deposition Modeling 


As material approaches the surface, deposition begins, and continues until 
all material is deposited (if evaporation occurs, some of the material will be 
left in the atmosphere to drift). Ground deposition is computed by assuming 
that the concentration of material around the mean may be taken as Gaussian 


exp[- G2) exp| - eee (44) 


20 90° 





ona2 


where the releaSed material is at- position (Y,Z) and | (y,z)” defines” its 
spatial effect. 


GAUSSIAN DEPOSITION 


The material is assumed to deposit entirely at the point of surface 
impact. Here, Eq. (44) is integrated across all z values to give 


2 

1 tt } 

Molane at 7 exp|- a 2 
¥210 20 





(45) 


CONTINUOUS DEPOSITION 


The material is assumed to deposit incrementally as it approaches the 
surface, and to continue depositing at its approach rate after impact. Here, 
Eqe (44) is integrated from negative infinity to the material location to give 


Z : 
= > =xp|= Ya) erfc eS (46) 


0° v20 


where erf@- is’ the complementary error function. The ground deposition is 
obtained by summing the incremental contributions to Molane as the 
integration proceeds. It may be seen that for material falling vertically 
toward the surface, the ground deposition pattern generated by Eq. (46) will 
be identical to the Gaussian deposition of Eq- (45). The use of Eq. (46) 
should lead to more realistic ground deposition distributions. 


CANOPY DEPOSITION 


When released material traverses a canopy, deposition on the canopy 
reduces the amount of material available for deposition on lower levels of the 
canopy or on the surface. The amount of material captured on the canopy 
vegetation will be directly related to the capture efficiency 8 and the 
distance traveled by the material through the canopy. Thus, at any time 
increment At , the deposited material may be found from 
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AC = BCA(z) Ad (47) 


wnere mad Pistithe distance’ traveled in At . The quantity AC reduces the 
remaining material C available for further deposition. 


Equation (47) contains A(z) p. the? plant area fraction, Thise is 
essentially a discrete function of height through the canopy, and is the ratio 
of leaf area (space “occupied" by the tree) divided by the surface area 
allocated to each typical tree. Reference 24 measured the leaf area for 
typical almond trees in the Chico orchard. Since the trees were spaced 8.2 
meters apart, the surface area for each tree becomes 68 m-* ‘The resulting 
plant area fraction, discrete at seven heights through the canopy, is shown in 
Figure 2-7. 


OBJECT DEPOSITION 


Equation (44) may also be interpreted to recover the deposition on 
collection devices (cards, cylinders or spheres) placed at specific locations 
in the wake flowfield. The total amount of deposition, as computed by 
Eq- (44), is reduced by the efficiency of the collection device, using the 
formulas and empirical data found in Ref. 25. The efficiency (as plotted in 
Figure 2-8 for the three devices) depends on two impaction parameters 


fe) pu 
im = ao (48) 
Pair slilie ae 


se eal 


fo en paki hawshae (49) 
Vair?p 


The significant target dimension D, is Ja typical cross-sectional distance 
sed to charaeterizesthe sizesof the™target (for.awcard, it is the smaller 
fength, for @ecylinder, it" ts? thessmaller tor wits diameter on length; for a 
sphere, it is its diameter). Figure 2-8 demonstrates that as the target 
dimension increases, K decreases from Eq. (48), 96 increases from Eq. (49), 
and the collection efficiency decreases (it is easier for a droplet to get out 
of the way of a big target that significantly disturbs the flowfield than a 
small target that doesn't). 


Material impaction also requires a knowledge of the orientation of the 
target relative to the solution coordinate system. For the purposes of this 
analysis, a normal vector must be supplied. This (x,y,z) vector defines the 
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Inertial impaction efficiencies for targets discussed in 
Ref. 25. The top curve in each graph gives 6 = 0 ; the 
middle curve gives ¢ = 100 ; and the bottom curve gives 


@ = 1000 . 


Figure 2-8. 


direction of a perpendicular line away from the target (in the case of a 
cylinder or a sphere this information isolates which half of the target is to 
be investigated). Some typical normal vectors are given in Figure 2-9. The 
effect of target orientation may be seen in Figure 2-10, where the deposition 
and target collection efficiency are plotted as functions of normal vector 
direction for a specific target location in Example Case 3. 





Overall Model Validation 


The first significant comparison of the AGDISP code predictions with data 
was reported in Ref. 26, where the influence of fully rolled-up tip vortices, 
propeller and crosswind in the mean and variance equations were favorably 
compared with a series of fixed-wing experiments. Later detailed comparisons 
(Refs. 27-30) demonstrate the favorable performance of the AGDISP code when 
compared with detailed field test data for both fixed-wing and helicopter 
experiments over a wide material size range. 


Zee 
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Figure 2-9. Typical target normal vectors given as (x,y,z) components. 
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Figure 2-10. The deposition and target collection efficiency on a sphere of 
diameter 0.075m located in the wake of the AgTruck (Example 
Case. 3) ah y= "—=10mm 5 Zo= 2m <. All normal wveetots are 
horizontal, with a direction of 0 degrees pointing in the -x 
direction (the flight direction of the aircraft); +90 degrees 
points in the +ty direction; -90 degrees points in the -y 
direction; and +180 degrees points in the +x direction (the wake 
of the aircraft). 


‘3. CODE OPERATION 


The AGDISP code consists of two programs: AGDISP for establishing the 
desired background fields and computing the material trajectories, and AGPLOT 
for plotting the resulting solutions. These programs are configured as 
interactive programs reading their needed input data from disk files and 
writing output to the terminal. Batch operation is discussed in Section 7. 
AGDISP reads an input file of data cards (valid cards are discussed in Section 
4), types information to the terminal as the run proceeds, generates printed 
output in a separate file, stores material trajectory information in a plot 
file for subsequent plotting, and reads WAKE plot file results from a separate 
data file (if invoked). AGPLOT reads the desired material plot file and 
interactively questions the user to supply necessary scaling data for the plot 
construction. 


The AGDISP code reads and processes the input data disk file, causing 
termination of the program if errors are discovered or needed cards are 
missing. Initialization follows: material locations and velocities, vortex 
strength, Betz roll-up, plant canopy, crosswind velocity profile, WAKE plot 
file. The material equations are then integrated repeatedly until one of 
three termination conditions occurs: 

l. maximum simulation time is reached; 


2 all of the released material has deposited on the surface; or 


36 all of the material has reached the evaporation cutoff diameter and 
are removed from the calculation (if invoked). 


Each integration step in turn includes the following: 
l. equation system solution of the equations of motion; 


Zs updating the background mean velocity and turbulence fields by the 
WAKE plot file, Betz roll-up, vortex motion, or plant canopy; and 


3. incremental file save of the results for subsequent plotting by 
AGPLOT. 


At the end of a run, the surface impact statistics are summarized and the 
deposition fraction (the mass fraction of released material that has reached 


the surface in the desired simulation time) is computed. 


AGDISP is run from and produces data into files assigned to specific unit 
numbers. These files are the following: 


3=1 


UNIT 


10 


Fl=25 


DESCRIPTION 


AGDISP input file containing the appropriately constructed data entry 
cards. 


Terminal input. 
Terminal output. 


CASEFILE default file set of AGDISP input files (accessed with input 
card 0005). 


Binary file constructed by AGDISP and plotted by AGPLOT. 
AGDISP printer listing file. 
WAKE binary plot file. 


Additional AGDISP binary plot files to be combined for composite 
deposition plots in AGPLOT. 


The current limits of the AGDISP code are the following: 


1. 


All of the idealized background mean velocity and turbulence fields 
are for neutral environments. The inputted crosswind velocity 
profile will be interpreted as locally neutral. 


The total number of nozzles that can be included in the simulation 
cannot exceed 60. The error message INCORRECT NUMBER OF DROPLETS 
will be invoked and AGDISP will terminate. 


The Betz roll-up cannot enter more than 100 discrete circulation 
values along the wing. The error message ERROR IN CIRCULATION DATA 
INPUT will be invoked and AGDISP will terminate. 


The Betz roll-up procedure cannot handle more than four discrete 
vortices rolling-up on each wing. In this case, the error message 
BETZ WILL ROLL-UP MORE THAN 4 VORTICES will be invoked and AGDISP 
will terminate. 


The plant canopy input data cannot permit more than 20 discrete 
entries to define the canopy vertical profile shape. The error 
message ERROR IN PLANT AREA INPUT will be invoked and AGDISP will 
terminate. 


The discrete crosswind velocity data cannot permit more than 20 


entries to define the velocity profile. The error message ERROR IN 
CROSSWIND VELOCITY INPUT will be invoked and AGDISP will terminate. 


Sey. 


the 


The wide body data cannot permit more than 20 entries to define the 
- fuselage cross-sectional area as a function of distance from the 


aircraft nose. The error message ERROR IN WIDE BODY AREA INPUT will 
be invoked and AGDISP will terminate. 


The input deck to AGDISP cannot exceed 200 cards. The message AGDISP 
CONFIGURED FOR 200 INPUT CARDS will be invoked and AGDISP will 
terminate. 


The complete alphabetized list of error messages generated in AGDISP are 
following: 


AGDISP CODE DOES NOT SUPPORT CARD 


the particular input card in question is not an allowable input card 
to AGDISP. 


AGDISP CONFIGURED FOR 200 INPUT CARDS 
a maximum of 200 input cards are permitted in the AGDISP input deck. 


BETZ WILL ROLL-UP MORE THAN 4 VORTICES 
a maximum of four vortices are permitted on a wing semispan. 


CARD ORDER INCONSISTENT AT CARD 
input cards are not in increasing order, or an error in input has 
lead AGDISP to expect a card that is not in the input deck. 


ERROR IN CIRCULATION DATA INPUT 
the Betz roll-up procedure must have between 3 and 100 discrete 
locations of circulation, with the last y location entered as the 
negative of its actual value. 


ERROR IN CROSSWIND VELOCITY INPUT 
the user-supplied crosswind velocity profile must have between 3 and 
20 discrete locations of velocity, with the last z height entered 
as the negative of its actual value. 


ERROR IN PLANT AREA INPUT 
the plant canopy must have between 3 and 20 discrete locations of 
plant area fraction, with the last z height entered as the negative 
of its actual value. 


ERROR IN WIDE BODY AREA INPUT 
the wide body cross-sectional area must have between 3 and 20 
discrete locations of area, with the last x distance entered as the 
negative of its actual value. 


FULL-PLANE CALCULATION REQUIRED 
Certain AGDISP options (terrain, propeller, crosswind, etc.) require 
that the solution be full-plane. If the user enters the half-plane 
flag on card 0010, AGDISP will terminate. 


INCORRECT NUMBER OF DROPLETS 
more than 60 nozzle locations are invoked in the simulation. 


INITIAL CONDITION MISMATCH 
Initial conditions were entered for a nozzle location that doesn't 
exist. 


INPUT DOES NOT FULLY INITIALIZE AGDISP 
input data cards are missing. 


INSUFFICLENT DATA BEFORE CARD 
the AGDISP run is not fully initialized; perhaps cards are missing 
from the input deck or a card does not have enough information on it. 


INVALID CASE NUMBER 
AGDISP could not locate the desired CASE number on card 0005 in the 
case file. 


PREMATURE END OF WAKE PLOT FILE REACHED 
the WAKE plot file ended without enough profile information. 


WAKE PLOT FILE EXTRAPOLATION 
Extrapolation of data in the WAKE plot file first occurs; this is a 
warning message. 


WAKE. PLOT FILE MESH SIZE 
the WAKE plot file must have between 2 and 16 mesh points in the y 
and .z ydvrections. 


The AGPLOT code interactively plots the AGDISP plot file contents on 
appropriate terminals or plotting devices. After scanning the file AGPLOT 
offers the current menu of plotting options available to the user (Section 5), 
who then determines the options to invoke. The list of error messages 
generated in AGPLOT are the following: 


MAXIMUM SCALE LESS THAN MINIMUM SCALE 
the user-supplied scale limits were inconsistent, and AGPLOT will 
Lequest. reentry of -the “scale data. 


MAXIMUM SCALE SIZE ADJUSTED 
the user-supplied scale limits were adjusted for consistency with the 
supplied scale delta. 


NO EQUIVALENT GAUSSIAN CONSTRUCTED 
there was no acceptable time at which an equivalent Gaussian could be 
constructed by AGPLOT. 


OPTION NOT AVAILABLE FOR THIS DATAFILE 


AGPLOT cannot invoke a_ requested option because the data is not 
available on the AGDISP plot file. 


3-4 


SCALE INCREMENT LESS THAN ZERO 


the user-supplied scale increment was inconsistent, and AGPLOT will 
request reentry of the scale data. 


TOO MANY SCALE DIVISIONS 


the user-supplied scale data forced more than ten scale divisions, 
and AGPLOT will request reentry of the scale data. 


TOTAL MASS FRACTION IS NOT 1.0 


the user-supplied mass fractions do not add to one for the multiple 
plot files examined. 


On the Data General, the following error messages are also possible: 


GKS_DRAW ERROR 
GKS OPEN ERROR 
GKS_Q WS _WIN ERROR 


GKS graphics errors on the Data General have occurred. These 
messages can only imply that something is significantly wrong with 
the program operation of AGPLOT. 
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'4. AGDISP INPUTS 


This section of the AGDISP User Manual details the input cards to the 
code. All data entry is in free format, with card data separated by commas or 
blank spaces. This convenience offers ease of formatting the data_but 


ee 


requires that every data card have all of its data values present, even if 
they are zero. Unless noted below as integer values, data is entered as real 


numbers (with decimal points). The MKS system is used throughout. 





All data cards begin with a four-digit identification in columns 1-4, with 
the rest of the data in free format after column 5. The order of the cards is 
important and must follow ever-increasing identification numbers. The AGDISP 
code has been programmed to verify this order. In addition, certain available 
options are inconsistent with each other. The AGDISP code has been programmed 
to trap these inconsistencies, with the message CARD ORDER INCONSISTENT AT 
CARD followed by the offending card when an error is detected. Unsupported 
cards are flagged with the message AGDISP CODE DOES NOT SUPPORT CARD, while 
missing data cards are flagged with INPUT DOES NOT FULLY INITIALIZE AGDISP. 
With free format there exists the chance that all necessary data does not 
appear on the appropriate card; in this case, the AGDISP code will run out of 
data cards before it rationalizes all of its data pointers. The error message 
will read INSUFFICIENT DATA BEFORE CARD following the card where data ended. 
The appearance of this error message may require a detailed examination of all 
data cards to determine the missing data. On the other hand, because all of 
the data is echoed on the terminal, error messages are traceable in a 
systematic manner. A careful check should be made the first time a new case 
is started, and at least until the user is confident of the program input 
requirements. Only by looking at what the program thinks is inputted will the 
user be able to verify that what was inputted was correct. In almost all 
cases, the code makes no check of input validity, either signs or 
magnitudes. The acceleration of gravity is built into the program as 
9.8 m/sec - Every AGDISP run requires the entry, at least, of cards 0010, 
0020, 0050 and 0060. Table 4-1 summarizes all input cards. 





Detailed description of the currently available data cards follows, 
including the special messages they invoke. 


Caution: all input cards are described below, but for a specific AGDISP 
run, only those cards that are needed to describe the run should be 
included. Default entry (cards 0010, 0020, 0050 and 0060) will be suf ficient 
to describe most runs. For example, if the run is made without evaporation, 
there is no reason to include card 0065. 


0000 CMNT 
Variable Type Description (Comment Card) 
CMNT CHAR Comment card, any number of which may be _ placed 


anywhere in the input deck. 


© 
= = 
on i= 
si BES 
Dub Bind cS ll seen cs 
By i. Sage eab is 
ona ESS ES 
mH DCOOmMd) & 
Card —_ Contents 
0000 Comment O 
0005 Casefile OF 
0010 Program @ . 
0011 Printer/Plotter Control Sth : 
0015 —- Terrain Slope Angle was . © 
0020 Aircraft Characteristics eo... 
0021 _ Biplane Characteristics OT 
0023. ~— Aircraft Weight d. 
0025 _—_— Betz Wing Load Distribution OD tere Ga 
0028 Neutral Crosswind «Qe. 
0029 ~— Discrete Crosswind Jenetans SQ ae 
0030 _—_— Helicopter aed jie 
0035  ——- Jet Engine See) 
0040 Propeller BD x 
0045 Powerplant Placement ~@ 
0050 Turbulence @ . yah 
0055 Canopy e . 
0056 Canopy Collection Efficiency oll Qit. 
0060 Nozzle Qieyrivoos 
0061 Discrete Nozzle Location EDAD A Sai 
0062 Nozzle Initial Condition RONG 13 0 feta por 
0065 Evaporation . @. 
0066 Evaporation Parameterization was Ee) a8 
0070 Apparent Surface Height Cae. 
0075 Circulation Decay POLLO) ee 
0080  Wide-Body Setup 1G) Oap 
008 1 Wide-Body Effect i) Oe 


@ Card must be present in input deck 
O Card may be present in input deck 
@@O8 = Card pairs, one of which must be present if invoked 


TABLE 4-1. AGDISP input card Summary, broken into convenient categories for 
typical options. The "Required" cards are added with other option 
cards to build the complete input deck. 
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0005 ICASE 


Variable Type Description (Casefile Card) 
ICASE INT Case to run (or modify) found in the default case 


input datafile. All of the input cards found on the 
case file for this case value will be added to the 
input runstream to build a complete input deck for 
AGDISP. Cards present in the input file will replace 
any cards of identical number found in the case file. 


The case file may contain any number of input files to 
AGDISP, each with its input information consistent 
with the discussion in this section. The only 
attribute that separates the case file decks is a 
first card to each deck containing the word CASE in 
columns 1-4 and the case number beyond column 5. With 
card 0005 present, AGDISP will search the default case 
input datafile to find a match of case number; if that 
match is made, all of the subsequent cards up to the 
next CASE card will be copied to the input runstream. 
If a match is not made, AGDISP will print INVALID CASE 
NUMBER and terminate. 


0010 TMAX LHFPL 

Variable Type Description (Program Card) 

TMAX REAL Maximum time in seconds for the simulation to run, a 
time which will be shortened if the material impacts 
the surface or is discarded at the evaporation cutoff 
before the maximum time is reached. 

LHF PL INT Simulation plane entry configuring either a half-plane 
(ve mO me SOlution 8(a .Vvalue of. °1)5 or. a full plane 
(- < y < +”) solution (a value of 2). The full plane 
solution should be invoked if crosswind, propeller or 
terrain exist. 

0011 NPTJ NPVX NPVL NPXX 

Variable Type Description (Printer/Plotter Control Card) 

NPTJ INT Modified default printing interval. Under normal 


operation AGDISP output is generated at appropriate 
intervals to the printer and the plot file. Entry of 
a number greater than 1 increases’ the _ printing 
interval (the plotting interval cannot be altered). 
Default value is 1; to turn off printing requires an 
entry of 0. 
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NPVX 


NPVL 


NPXX 


OU1L5 


Variable 


TA 


TA 


INT 


INT 


INT 


Type 


REAL 


Invokes plot file saving of the time histories of the 
driving trajectories, and controls the interval of 
printing these time histories. The driving 
trajectories are those associated with the centers of 
vortices, the center of the propeller and the center 
of the helicopter rotor disk (whenever any are present 
in a particular run). Any value for this entry 
invokes additional plot file output, to permit 
subsequent plotting of the driving trajectories by 
AGPLOT. A number greater than | increases’ the 
printing interval. Default value is 0. 


Controls whether the vertical velocities are saved on 
the plot file for subsequent plotting by AGPLOT. 
Default value is 0 (no saving of vertical velocities); 
any nonzero value entered on this card will force 
additional plot file output. Vertical velocities must 
be saved if the AGPLOT option for the equivalent 
Gaussian distribution is to be invoked. 


Controls whether the mean velocities are saved on the 
plot file for subsequent plotting by AGPLOT. Default 
value is 0 (no saving of mean velocities); any nonzero 
value entered on this card will force additional plot 
file output. Mean velocities must be saved if the 
AGPLOT option for deposition on objects is to be 
invoked. 


In all cases, additional writes to the plot file will 
increase its size. 


Description (Terrain Slope Angle Card) 


Local ground slope in degrees. A positive terrain 
angle raises the right side of the surface as viewed 
from behind the aircraft, with the origin remaining 
along the initial aircraft vertical centerline. The 
terrain is characterized by a locally straight surface 
so that all of the simplified flow field options in 
AGDISP remain available. A full plane solution is 
necessary. Canopy and crosswind effects remain 
parallel to the tilted surface. 
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0020 


Variable 


LMVEL 


wG 


DIST 


UO 


0021 


Variable 


DZBP 


PSBP 


PGBP 


LMVEL S_ DIST 


Type 


INT 


REAL 


REAL 


REAL 


DZBP PSBP 


Type 


REAL 


REAL 


REAL 


UO 
Description (Aircraft Characteristics Card) 


Mean velocity flag taking one of seven values: 


5 denotes ground sprayer; 
4 denotes helicopter entry; 
3 


denotes an elliptically loaded, fully rolled-up 
tip vortex; 
2 denotes a rectangularly loaded, fully rolled-up 


tip vortex; 
1 denotes a triangularly loaded, fully rolled-up tip 
vortex; 
0 denotes Betz roll-up from a given circulation, and 
-]1 denotes WAKE plot file entry (explained under card 
0050). 


An entry of 4 requires a 0030 card; entry of 3, 2 or 1 
requires a 0023 card; and an entry of O requires 
multiple entries of 0025 cards. 


Semispan of the aircraft in meters 
radius for a helicopter). 


Calso. the, rotor 


Nominal height of the aircraft wing about the surface 
in meters (this distance is the assumed nominal nozzle 
release height). For fully rolled-up tip vortices, 
this height is the initial Z coordinate of the 
vortex centerline. For Betz roll-up data this height 
is the Z coordinate of the initial vortex sheet 
(assumed horizontal). 


Flight speed of the aircraft in m/sec. 


PGBP 


Description (Biplane Characteristics Card) 
Vertical distance in meters from the main wing 
location specified on card 0020 to the biplane wing. 


Semispan of the biplane wing entered as a fraction of 
the semispan of the main wing (if the wings are equal 
in length this entry would be 1.0). 


Weight carried by the biplane wing entered as a 
fraction of the weight carried by the main wing (for 
equal weights this entry would be 1.0). 


0023 


Variable 


WT 


0025 


Variable 


Y 


G 


0028 


Variable 


U 


G 


ZO 


CA 


WT 


Type 


REAL 


eG 


Type 


REAL 


REAL 


Of 4 780) (GA 
Type 
REAL 
REAL 
REAL 


REAL 


Description (Aircraft Weight Card) 


Weight of the fixed-wing aircraft in N , for a fully 


rolled-up tip vortex pair. 


Description (Betz Wing Load Distribution Cards) 


Position in meters measured from the root 


monotonically towards the tip. 


wing 


Circulation value in m2/sec at .this. location. 


The final card is signaled by entering the last y 
position as the negative of its actual value. This 
card then forces initialization of the Betz roll-up 
procedure. This initialization (and the ensuing roll- 
up) invoke extra printer output summarizing the roll- 
up process. 


As the Betz roll-up continues, the material being 
tracked will be influenced by ever-increasing strength 
vortices, whose positions and strengths will approach 
those of fully rolled-up vortices. The ACDISP code 
includes the effect of the unrolled-up sheets on the 
motion of the released material. 


Description (Neutral Crosswind Card) 


Mean wind velocity in the (x,y) plane in m/sec. 
Altitude of the mean wind velocity in meters. 
Surface roughness height 25 in meters. 

Direction angle in degrees from which the wind is 
blowing, i.e., for 0° the wind is a head wind; for 
180°. the wind is a tail wind; for 90° the wind is a 


crosswind from right to left; and for -90° (or 2709 
the wind is a crosswind from left to right. 
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0029 Z SUIS 
Variable Type 
6 REAL 
U REAL 
V REAL 
0030 WT BDOT 
Variable Type 
WT REAL 
BDOT REAL 
0035 THT RPRP 
Variable Type 
THE REAL 
RPRP REAL 
DZ REAL 


DZ 


Description (Discrete Crosswind Cards) 


Position in meters beginning at the surface and 
increasing monotonically to the top of the 
computational domain. 


Axial velocity value in m/sec, positive pointing 
downstream (headwind). 


Crosswind velocity value in m/sec positive pointing in 
the positive y direction. The square root of the 
sum of the squares of U and V_ recovers the total 
velocity magnitude at height z. 


The final card is signaled by entering the last Zz 
position as the negative of its actual value. This 
card entry then initializes a calculation of the 
turbulence level throughout the velocity profile. If 
a canopy is present, it is assumed that the entries on 
cards 0029 have been modified by the user of AGDISP to 
reflect the influence of the canopy on the crosswind 
velocity profile and_ the subsequent turbulence 
calculation. 


Description (Helicopter Card) 


Weight of the helicopter in N. 
Blade rotation rate in rpm. 
The helicopter is idealized as a dynamic transition 


from a rotor downwash field to a rectangularly loaded, 
fully rolled-up vortex pair. 


XPRP 


Description (Jet Engine Card) 


Thrust of the engine in N. 
Exit radius of the engine in meters. 


Incremental distance in meters of the engine 
centerline above or below the nominal release height 
given on card 0020 (the jet engine is assumed to be at 
the airplane centerline, y = 0)- The placement of 


XPRP 


0040 CD 


Variable 
CD 

AS 

ETA 

TDOT 
RPRP 


DZ 


XPRP 


0045 NPRP 


Variable 


NPRP 


TV 


REAL 


the jet engine is superceded by the presence of card 
0045. 


Axial distance from the trailing edge of the wing to 
the exit plane of the jet engine in meters, where 
positive is downstream. 


AS z*ETA™ “TDOT “RPRP: -DZePXPRP 


Type 
REAL 
REAL 
REAL 
REAL 
REAL 


REAL 


REAL 


Description (Propeller Card) 


Drag coefficient of the aircraft. 


Airplane planform area in m?. 


Propeller efficiency. 
Shaft rpm. 
Propeller blade radius in meters. 


Incremental distance in meters of the shaft centerline 
above or below the nominal release height given on 
card 0020 (the propeller is assumed to be at the 
airplane centerline, y = Q). The placement of the 
propeller is superceded by the presence of card 0045. 


Axial distance from the trailing edge of the wing to 
the propeller blade plane in meters, where positive is 
downstream. 


(TV(N), N = 1 TO NPRP) 


Type 


INT 


REAL 


Description 


(Powerplant Placement Card) 


Number of engines (1 to 4, although | is the default 
on cards 0035 and 0040). 


An equal number of dimensions following the pattern 


for one engine: 
incremental distance (meters) at the aircraft 
centerline above or below the nominal height. 


for two engines: 
positive horizontal position (meters) of the engine 
from the aircraft centerline; and 
incremental distance (meters) of the engine above 
or below the nominal height (the other engine is 
Symmetrically positioned); 
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for three engines: 
incremental distance (meters) at the aircraft 
centerline above or below the nominal height for 
the engine at the aircraft centerline; 
positive horizontal position (meters) of the second 
engine from the aircraft centerline; and 
incremental distance (meters) of the second engine 
above or below the nominal height (the third engine 
is symmetrically positioned); 


for four engines: 

positive horizontal position (meters) of the first 
engine from the aircraft centerline; 

incremental distance (meters) of the first engine 
above or below the nominal height (the second 
engine is symmetrically positioned); 

positive horizontal position (meters) of the third 
engine from the aircraft centerline; and 
incremental distance (meters) of the third engine 
above or below the nominal height (the fourth 
engine is symmetrically positioned). 


0050 LQQSE QQMX 
Variable Type Description (Turbulence Card) 
LQQSE INT Turbulence flag carrying one of the following values: 
1 invokes superequilibrium turbulence; 
0 assumes a fixed value of turbulence given on this 
card; and 
-l1 specified the attached WAKE plot file (invoked 
with a -l entry on card 0020). 
QQMX REAL Maximum value of the background turbulence q? in 


m /sec* (this value is incremented by the presence of 
a crosswind). The turbulence level generated by a 
crosswind flowfield is generally sufficient to 
represent the ambient turbulence level (the input on 
card 0050 would then be O.O for q2). However, in 
some cases, it may be necessary to augment this value 
for specific local atmospheric conditions. A short 
guide to the selection of appropriate background 
turbulence levels is offered in Section 8. 


A -1 on card 0020 invokes the WAKE plot file and must 
be accompanied by a -1 on this card. The WAKE plot 
file contains the crossplane velocities V and W , 
and the turbulence q“~ . 


0055 


Variable 


Z 


Z 


AA 


Type 


REAL 


REAL 


The sequential binary WAKE plot file is constructed as 
follows: 


record 1: the number of y (or horizontal) mesh 
points and z (or vertical) mesh points in the plot 
file (two integers). The code restricts both of these 
entries to 16 or less; larger walues invoke WAKE PLOT 
FILE MESH SIZE; 


record 2: the y mesh values in meters; 
record 3: the z mesh values in meters; 
record 4 and following: the profile data follows, 


with each time slot repeating the same pattern. It 
begins with a single record which is the time saved in 
seconds. The data for each variable on the file 
(specified in the order V,W,q2 ) follows, by giving 
all, +of,) the y values of V at the first Zz 
position, then, all) ofthe, jy values of V_ at the 
second z position, on to the last z _ position; then 
onto sadli ofe them ys tvalues ‘of W at the first ag 
position, etc<; suntil-~ ‘the «values. of “all “of wtne 
variables at all of the y and 2z _ positions have 
been given. The next time value follows, repeated to 
the end of the plot file. Velocities are in m/sec 
and turbulence in m/sec“ . The warning message WAKE 
PLOT EILE. EXTRAPOLATION,.is, ,output they first. tine 
spatial extrapolation must be used by the AGDISP code 
during interpolation for the variables contained on 
the WAKE plot file. For times beyond the entries in 
the WAKE plot file, the spatial profiles nearest in 
time to the time being solved in AGDISP will be used. 





Before integrating the equations, AGDISP reads the 
entire WAKE plot file through to its end. If the file 
is constructed incorrectly, the error message 
PREMATURE END OF WAKE PLOT FILE REACHED is invoked, 
and AGDISP terminates. 


Description (Canopy Cards) 


Position in meters beginning at the surface and 
increasing monotonically to the top of the canopy. 


Plant area fraction corresponding to the 2z_ position. 


The final card 0055 is signaled by entering the last 
z position as the negative of its actual value. This 


card entry then initializes the canopy calculation, 
computing the displacement thickness of the canopy and 
forcing modification to the crosswind velocity and 
turbulence within the canopy. During the AGDISP run, 
the trajectories of the aircraft vortices will be 
altered upon entering the canopy, and the position of 
the helicopter dividing streamline will be changed. 





0056 CEFF 

Variable Type Description (Canopy Collection Efficiency Card) 

CEFF REAL Canopy capture efficiency 8 . Default value is 1.0. 

U060 LPART LZERO DZ XO DIAM ODENF 

Variable Type Description (Nozzle Card) 

LPART INT Number of nozzles in the half-plane (y> 0). 

LZERO INT Centerline »nozzle flags’>equal®’to 1° if a nozzle is 

positioned at the centerline of the aircraft (y = 0); 
otherwise it is equal to 0. 
A half-plane solution will track total released 
material equal to the sum of these two integers; a 
full-plane solution will track total released material 
equal to the sum of twice the first integer plus the 
second integer. Additionally, if these two integers 
are entered positive, AGDISP will position the nozzles 
uniformly along the wing (for a first integer entry of 
1, the nozzle will be positioned at 1/2 the semispan; 
for a first integer entry of 2, the nozzles will be 
positioned at 1/3 and 2/3 the semispan, etc.). If 
these two integers are entered negative, AGDISP will 
expect sufficient 0061 cards in the input deck to 
initialize all of the nozzles. 

DZ REAL Vertical position in meters off-setting the material 
release point from the height of the wing given on 
card 0020. 

XO REAL Nominal axial position of the nozzles relative to the 
trailing edge of the wing (or the shaft centerline of 
the helicopter) in meters, positive measured 
downstream. 

DIAM REAL Diameter in microns of the released material (only one 


size may be released during any one AGDISP run). 


DENF 


006) 


Variable 


4 


Z 


0062 


Variable 


N 


XS 


VS 


REAL 


Type 


REAL 


REAL 


Nos Uy VaQw axs 


Type 


INT 


REAL 


REAL 
REAL 


REAL 


REAL 


Specific gravity of the released material. 


Description (Discrete Nozzle Location Cards) 
Y value in meters along the spray boom. 


Vertical position in meters off-setting the release 
point. from the height of the wing given on card 
0020. The vertical position off-set on card 0060 is 
disregarded. Discrete nozzle locations allow for fine 
tuning of nozzle locations, and are necessary if 
nozzles are not uniformly distributed along the wing 
or if the spray boom is not parallel to the surface. 


There must be enough 0061 cards to satisfy the 





requirements of the 0060 card. If a full-plane 
solution is invoked all nozzle locations must be 
specified. If a nozzle is located along the 


centerline in a half-plane simulation, that nozzle 
must be the last card 0061 for purposes of computing 
the ground deposition correctly. 

VS 


Description (Nozzle Initial Condition Cards) 


Nozzle number, recorded from the left most nozzle 
position. 


Initial axial velocity in m/sec at nozzle N 
positive measured downstream. 


Initial horizontal velocity in m/sec at nozzle N. 
Initial vertical velocity in m/sec at nozzle N. 


Initial spatial variance of the material path in m2 at 
nozzle N. 


Initial velocity variance of the material in m/sec 
at nozzle N. 


All initial conditions are set to zero unless modified 
by card 0062 for each nozzle. Card 0062 is ideal for 
configurations with few nozzles, such as ground 
sprayers. 
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0065 


Variable 


DTEMP 


DCUT 


0066 


Variable 


DEA 


DEB 


DEC 


0070 


Vartable 


4 


DTEMP DCUT 


DEA 


Z 


Type 


REAL 


REAL 


DEB DEC 


Type 


REAL 


REAL 


REAL 


Type 


REAL 


Description (Evaporation Card) 
Wet bulb temperature difference in deg C. 


Material diameter in microns below which evaporation 
ceases (the cut-off diameter). 


If the wet bulb temperature difference entered on this 
card is zero, AGDISP assumes that evaporation will be 
parameterized, and expects card 0066. LES DCUT as 
negative, material is removed from the simulation when 
it evaporates down to the absolute value of DCUT. 


Description (Evaporation Parameterization Card) 


Coefficient A in microns in the material diameter 
equation 


Diameter = A+ Bt + cr2 (50) 


Coefficient B in microns/sec. 


Coefficient +CeLin microns/sec@. 


These coefficients are compatible with the evaporation 
model in FSCBG (Ref. 31), and produce consistent 
evaporation results between AGDISP and FSCBG. The 
coefficient A should be identical to the material 
diameter entered on card 0060. 


Description (Apparent Surface Height Card) 


Apparent height above the surface (in meters) at which 
the material deposits. The default condition requires 
impact of material at the surface (with an apparent 
height equal to zero). For a nonzero entry, when 
material intersects the apparent height, it is removed 
from the computation, just as though it had impacted 
the surface. A subsequent ground deposition plot with 
AGPLOT will recover the deposition of the material at 
the apparent height entered on this card. 


0075 GDK 


Variable Type 
GDK REAL 
0080 XBOD DZ 
Variable Type 
XBOD . REAL 
DZ REAL 
0081 XvSS 
Variable Type 
X REAL 
mohens REAL 


Description (Circulation Decay Card) 


Factor (nondimensional) that forces a time reduction 
in the circulation strength of the vortices. Default 
value is zero (no decay). If *this entry is negative, 
its absolute value is taken as the circulation decay 
constant multiplied by the ambient turbulence level. 
Reference 14 suggests that for small material sizes, 
an appropriate entry on card 0075 is -0.56 m/sec. 


Description (Wide-Body Setup Card) 


Position in meters of the front of the fuselage from 
the trailing edge of the wing (or the shaft centerline 
of the helicopter), positive measured downstream (by 
convention, then, this number must be negative). 


Vertical distance in meters off-setting the nose 
centerline of the fuselage from the nominal height of 
the wing (or rotor plane) given on card 0020. The 
fuselage is assumed to be axisymmetric. If card 0080 
is present, cards 0081 must also be present to define 
the fuselage body shape. 


Description (Wide-Body Effect Cards) 


Position in meters, measured positive from the nose of 
the fuselage to the tail, increasing monotonically. 
Cross-sectional area of the in m2, 
corresponding to the x _ position. 


fuselage 


The first 0081 card must contain 0.0 for x and 0.0 
for the area. The final card is signaled by entering 
the last x position as the negative of the actual 
fuselage body length, with the corresponding area as 
0.0. This card entry initializes the wide body 
effects of a source-sink velocity profile near the 
fuselage and an axisymmetric turbulent wake decaying 
behind the fuselage. 


22.22 AGPLOT  INRUTS 


The AGPLOT program may be invoked to plot the resulting trajectories and 
deposition pattern from an AGDISP run or multiple AGDISP runs. The code asks 
questions of the user in the interactive environment, and with this 
information processes the data and plots the results. 


At the beginning of AGPLOT, the data file is read, the input data deck 
used to generate the file is written to the screen (to serve as verification 
of the proper run to plot), and the file is scanned for maximum values and 
available options. With this information in hand, AGPLOT will display its 
current menu driver: 


EXIT AGPLOT 

MEAN DROPLET TRAJECTORIES 

MEAN + STANDARD DEVIATION TRAJECTORIES 
VORTICES/HELICOPTER/ENGINE CENTROIDS 
GAUSSIAN GROUND DEPOSITION 

CONTINUOUS GROUND DEPOSITION 

CANOPY DEPOSITION 

TOTAL CANOPY DEPOSITION 

EQUIVALENT GAUSSIAN DISTRIBUTION 
CROSSWIND VELOCITY PROFILE 

10 CANOPY PLANT AREA FRACTION PROFILE 

11 DROPLET DIAMETER TIME HISTORY 

12 DROPLET CANOPY PARAMETER TIME HISTORY 
137 DROPLET AXIAL VELOCITY TIME HISTORY 

14 DROPLET HORIZONTAL VELOCITY TIME HISTORY 
15 ~ DROPLET VERTICAL VELOCITY TIME HISTORY 
16 OBJECT DEPOSITION TIME HISTORY 

17? DRTFT PRACTIONSIIME SHISTORY 


OMANAD UF WN OO 


followed by the request ENTER OPTION TO RUN. Only those options available 
from the data in the plot file will be displayed. Each of the options is 
described as follows: 


0 is entered to exit from AGPLOT. Since several menu options may 
be invoked for one AGDISP run, OQ is the entry used to complete 
the plotting for a particular plot file. 


1 is entered to plot the mean trajectories. 


2 is entered to plot the mean trajectories (as solid curves) and 
the standard deviation of the trajectories (as dashed curves). 


3 is entered to plot the trajectories of the vortex centroids (as 
solid curves) and the center(s) of the propeller, jet engine or 
helicopter (as a dashed curve), if applicable. This option is 
available only with an appropriate nonzero entry on card 0011. 


is entered to compute and plot the composite ground deposition 
pattern for the current plot file and for additional plot files 
(up to a total of 16 plot files) established by the user and 


accessed with unit numbers of 11 and higher. The user must 
decide the appropriate mass fractions for each of the runs to be 
combined. Since. the) currents,plot. fileyis» already jinw place, 


AGPLOT will request ENTER CURRENT FILE MASS FRACTION. After this 
AGPLOT will open unit 11 and request ENTER MASS FRACTION for the 
next plot file. Further plot files will be expected sequentially 
on ,units,12,5 13 ,el4,.ete. .eAs totalamass, fraction .of. unity wide 
terminate the requests. If the total mass fraction is not close 
to unity, the warning message TOTAL MASS FRACTION IS NOT 1.0 will 
appear. AGPLOT will represent all of the impacted particles as 
Gaussian distributions, following Eq. (45). Six different 
deposition scales are available in AGPLOT. These scales require 
the mass flow rate per nozzle, and so AGPLUT requests ENTER 
CALLONS/MINUTE/NOZZLE to properly compute the deposition. These 
six .. scales ware: NORMALIZED, LITERS/HECTARE, OUNCES/ ACRE, 
GALLONS/ACRE, MG/M**2 and DROPS/CM**2. AGPLOT requests ENTER 
SCALE TO USE. An integral under the deposition curve will 
reflect the number and mass of material that have impacted the 
surface. This number appears in a corner of the plot. 


is entered to compute and plot the composite ground deposition 
pattern as for option 4 but with continuous deposition following 
Eqs (46). 


is entered to compute and plot contours of the composite canopy 
deposition pattern following incremental Gaussian deposition 
through the canopy. Up to a total of 16 plots files may be 
combined in one contour plot, with the procedure similar to 
option 4. Here, however, the droplet canopy parameter defined in 
Eq- (47) is summed to accumulate deposition in specified layers 
through the canopy. The user is then quizzed to give specific 
contour values to be plotted. Option 6 must be exercised to 
generate the profile distribution needed for option 7. Contour 
values appear in a corner of the plot. 


is entered to plot the total canopy deposition. 


is entered to compute and display the equivalent Gaussian 
distribution. Since the position and standard deviation of all 
material in the simulation is computed by AGDISP as a function of 
time after release, the equivalent mean position and standard 
deviation may be determined by appropriate summation § and 
integration at each time over all material in the simulation. A 
measure of the compatibility of the equivalent Gaussian 
distribution is made by computing a figure of merit ranging from 
O- stouel.. When the equivalent Gaussian no where represents the 
multiple-particle distribution, the figure of merit equals 
zero. When the equivalent Gaussian is every where identical to 
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the multiple distribution, the figure of merit is unity. 
Equivalent Gaussian calculations are made for every time on the 
pPlot® file sw with Wall) vtdmes! | givaniva monotonically increasing 
sequence number. The critical sequence numbers, the time 
associated with each of them, and their respective figure of 
merit and terminal velocity criterion are typed to the screen 
whenever: the first sequence number is reached; the last sequence 
number is reached; the figure of merit reaches a local minimun; 
the figure of merit reaches a local maximum; or material first 
comes within a standard deviation of the ground. 


The terminal velocity for the material size examined is computed 
assuming no evaporative effects and a_ specific gravity of 
unity. It is displayed with the message TERMINAL VELOCITY as 
reference to the user. The terminal velocity criterion is 
computed as the average over all material of their vertical 
velocities relative to their terminal velocity (if all material 
were moving at their terminal velocity, the WCRIT value would be 
one). Option 8 requires that the particle vertical velocity be 
saved on the plot file by the appropriate input on card 0011 in 
AGDISP. Material whose vertical velocities are not within fifty 
percent of their terminal velocity will not contribute to the 
figure of merit calculations at that time step. 


After processing the entire plot file and displaying the above 
summary information, the plot file is reread to locate the 
position of maximum figure of merit. At thissipoint, «the 
equivalent Gaussian data is displayed on the screen. If the 
equivalent Gaussian vertical standard deviation contacts’ the 
surface, the warning message GROUND ENCOUNTER BY EQUIVALENT 
GAUSSIAN will be displayed on the screen. Material that does not 
contribute to the equivalent Gaussian distribution (because it is 
trapped by vortex flow in the aircraft wake, or because it has 
come within a standard deviation of the surface) will be noted by 
the **TRAPPED** comment. If none of the material in the 
Simulation qualify for contribution to the equivalent Gaussian, 
the message NO EQUIVALENT GAUSSIAN CONSTRUCTED will be displayed, 
and the subsequent plot will be skipped. 


The resulting plot represents each released material distribution 
as dashed curves and the single equivalent Gaussian distribution 
as solid curves. Around each center point of the material 
contributing to the equivalent Gaussian, normalized contour lines 
are plotted at one-fourth and one times the standard deviation. 
Material not contributing to the equivalent Gaussian are plotted 
with much smaller ellipses to aid in identifying them. The 
plotted figure of merit appears in a corner of the plot. 


is entered to plot the crosswind velocity profile (if available 
with a 0028 card or 0029 cards). 


is entered to plot the canopy plant area fraction profile (if 
available with 0055 cards). 
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is entered to plot the time history of a droplet diameter (if 
evaporation is occurring). If there are more than one droplet in 
the simulation, AGPLOT will request ENTER DROPLET TO PROCESS. 


is entered to plot the time history of the canopy deposition of a 
droplet (if canopy cards 0055 are present). 


is entered to plot the time history of the axial U velocity of 
a droplet (if the mean velocity data is stored on the plot file 
with card 0011 in AGDISP). 


is entered to plot the time history of the horizontal V 
velocity of a droplet (if the mean velocity data is stored on the 
plot file with card 0011 in AGDISP). 


is entered to plot the time history of the vertical W velocity 
of a droplet (if the mean velocity data is stored on the plot 
file with card 0011 in AGDISP). 


is entered to compute the time history of the deposition on a 
specified collector geometry at a specified position in the 
AGDISP solution field (the mean velocity data must be stored on 
the plot file with card 0011 in AGDISP). 


The collector must first be placed in the solution field by 
responding to the request ENTER (Y,Z) COLLECTOR LOCATION with 
the y and z values in meters. Next, AGPLOT asks for the type 
of collector; currently, the options are card, cylinder or 
sphere. Then AGPLOT requests ENTER SIGNIFICANT TARGET 
DIMENSION to give a typical length in meters of the collector. 
Finally, AGPLOT needs to know the orientation of the collector 
relative) toy.the Cxpy' zi) coordinate system, asking ENTER 
(X,Y,Z) NORMAL VECTOR. Here all three components must be 
entered; AGDISP will renormalize into a unit vector. A preferred 
normal direction must be supplied, even for collectors that do 
not appear to have planar geometry. If the collector were placed 
parallel to the surface, the normal vector would be 0,05 PS eetr 
the collector were vertical and deposition on the left side were 
needed, the normal vector would be (0,-1,0). 


As with other deposition options, up to 16 plot files may be 
combined to compute the collector deposition. The total 
collection efficiency appears in a corner of the plot. 


is entered to compute the time history of the drift of all of the 
droplets in the simulation (up? tomlemiplotuefiless* collected 
together). 


Where appropriate, the following additional requests will be made: 


FULL PLANE PLOT is asked whenever the AGDISP run is a half-plane run. 
Response in Y or N for yes or no, respectively. 
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AUTOSCALE AXES will compute appropriate plotting scales if invoked. If N 
is the response, AGPLOT will ask for Y SCALE MIN, MAX, INCR and Z SCALE 
MIN, MAX, INCR. AGPLOT has a scale checking routine and will invoke a 
warning message whenever the scale delta is not an integer fraction of the 
overall scale size. Additionally, three errors will be trapped: MAX SCALE 
LESS THAN MIN SCALE; SCALE INCR LESS THAN ZERO; and TOO MANY SCALE 
DIVISIONS whenever more than ten scale divisions are needed. 


AUTOSCALE GRID will compute appropriate deposition intervals, if invoked, 


for options 5 and 6. If N is the response, AGPLOT will ask for Y SCALE 
MIN, MAX. 


AUTOSCALE CONTOURS will contour plot appropriate deposition levels through 
the canopy in option 6. LE N is the response, AGPLOT will ask for 
contour levels (and will plot them) until a negative deposition value is 
entered. 


TAG LOCATION AT TIME INCREMENTS is asked for mean trajectories, to tag 
released material positions on the plot as a function of time. The size 
of the tag is scaled to the size of the material standard deviation. 5 
the response is Y, AGPLOT will ask ENTER TAG TIME INCREMENT to enter the 
incremental seconds between tags. 


ENTER PLOT TITLE is asked to enter the title of the current plot. 
On the Data General the following additional messages are possible: 


ENTER MASTER FILENAME 
The filename containing the data to be manipulated (8 characters or 
less without extension) is entered here. 


ENTER METAFILE.EXT FOR THIS PLOT 
If metafiles are open (from a previous query), then each plot must be 
assigned a name and extension for storing on the DG system in the 
user directory. A metafile extension of GKM is recommended for entry 
into CEO. 


ENTER PLOT SCALING FACTOR 
Typically, all plots produced by AGPLOT will fill to the size 
available in the device. By entering a number here less than one, 
the plot will be scaled accordingly. 


ENTER QUE NAME CODE 
A hardcopy unit must be further specified by its queue name within 
the DG environment. 


ENTER WORKSTATION TYPE 
Response to this query will determine the plotting device to be 
used. If the first character of the response is not "D° or ‘d", 
AGPLOT assumes the device is a hardcopy unit and will request its 
QUE NAME. 
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OPEN OUTPUT METAFILES 
While working with any graphics device, the resulting plots may also 


be spooled to a file for later manipulation. This option permits 
that operation. 


PLOT COMPLETED: PRESS *RETURN* TO CONTINUE 


The DG system will blank the screen after completing a plot unless a 
pause is programmed into AGPLOT. 


Additionally, on the DG system, the pressing of any function key at a query 
line will return the user to the menu driver. 


6. TEST CASES 


Four test cases are included in the User Manual to illustrate the variety 
of data entry available with AGDISP. These examples are included to demon- 
strate the structure of the input card deck and the data entry requirements of 
each type of data card and option. They are not meant to suggest that certain 
options go with certain other options, nor are all options covered. 


The four examples illustrate the following features: 


l. 


The first test case examines the wake behind a F-15 fighter flying at 
Mach 0.5. Because of the high speed, the Betz roll-up is invoked, 
with a wing load distribution obtained from vortex-lattice theory. 
Figure 6-] displays the appropriate input deck. Figure 6-2 shows the 
mean material trajectories, and Figure 6-3 gives the Gaussian 
deposition for a release of 1 gallon/minute/nozzle. 


The second test .case examines the wake behind a wide body C-130 
transport plane. Figure 6-4 displays the appropriate input deck. 
Figure 6-5 shows the mean material trajectories for the assumed 
nozzle release points (the four propellers skew the solution); Figure 
6-6 shows the paths of the tip vortices and the four propellers, and 
Figure 6-7 shows the continuous ground deposition for a release of 1] 
gallon/minute/nozzle. 


The third test case examines the wake behind a Cessna 188 AgTruck in 
forward flight (Mission test run 14, droplet size 171 microns, Ref. 
30). Figure 6-8 displays the input deck while Figure 6-9 shows the 
mean material trajectories and Figure 6-10 gives the Gaussian ground 
deposition for a release of 0.1] gallon/minute/nozzle. Figure 6-11 
illustrates the equivalent Gaussian distribution, Figure 6-12 shows 
the time history of the evaporating diameter of one droplet in the 
Simulation, Figures 6-13, 6-14 and 6-15 give the droplet velocity 
time histories, and Figure 6-16 shows the deposition on a spherical 
collector. 


The fourth test case examines a Hiller 12E helicopter (Chico test run 
B-2, Ref. 27). Figure 6-17 displays the input file, while Figure 
6-18 shows the mean material trajectories. Figure 6-19 displays 
contours of canopy deposition, while Figure 6-20 summarizes the total 
deposition in the canopy, on the same scales as the continuous ground 
deposition, Figure 6-21. Figure 6-22 illustrates the crosswind 
velocity profile, Figure 6-23 illustrates the canopy plant area 
fraction, Figure 6-24 gives the canopy deposition fraction for a 
droplet in the simulation, and Figure 6-25 summarizes the drift 
fraction for the simulation (since there is no evaporation, the final 
drift fraction will reach a value of zero). 
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Figure 6-4. AGDISP input file for Example Case 2. 
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Figure 6-8. AGDISP input file for Example Case 3. 
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0000 CHICO B-2: HILLER 12E 
0010 150.0 2 

0011 100 0 0 0 
0020 4 5.4 14.20 11.09 
0028 1.2 16.0 0.1 67.0 
0030 10325.0 400.0 
0050 0 0.0 

0055 0.0 0.044 
0055 0.6 0.044 
0055 1.8 0.147 
0055 3.1 0.221 
0055 4.3 0.191 
0055 5.5 0.150 
0055 6.7 0.074 
0055 -7.3 0.0 
0056 0.2 

0060 -14 0 0.0 0.0 251.0 1.0 
0061 -4.12 -2.0 
0061 -3.81 -2.0 
0061 -3.51 -2.0 
0061 -3.20 -2.0 
0061 -2.90 -2.0 
0061 -2.59 -2.0 
0061 -2.29 -2.0 
0061 -1.98 -2.0 
0061 -1.68 -2.0 
0061 -1.37 -2.0 
0061 -1.07 -2.0 
0061 -0.76 -2.0 
0061 -0.46 -2.0 
0061 -0.15 -2.0 
0061 0.15 -2.0 
0061 0.46 -2.0 
0061 0.76 -2.0 
0061 1.07 -2.0 
0061 1.37 -2.0 
0061 1.68 -2.0 
0061 1.98 -2.0 
0061 2.29 -2.0 
0061 2.59 -2.0 
0061 2.90 -2.0 
0061 3.20 -2.0 
0061 3.51 -2.0 
0061 3.81 -2.0 
0061 4.12 -2.0 


Figure 6-17. AGDISP input file for Example Case 4. 
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7. ‘USING AGDISP AND AGPLOT 


AGDISP and AGPLOT are operational on several computer systems. For 
completeness, the Univac 1108 system at Ft. Collins will also be discussed, 
although AGDISP Mod 4.0 is its highest level of operation. 


Operation on the Data General MV-15000 


AGDISP and AGPLOT are run on the Data General Forest Service network by 
the commands 


XEQ AGDISP 
XEQ AGPLOT 


with resolution of the file names and extensions performed within the 
programs. 


The plotting package GKS is invoked to plot AGPLOT results onto graphics 
output devices. 


Operation on an IBM PC/XT/AT 


AGDISP on a PC is detailed in Ref. 32, with only a summary given here. 
The DOS operating system does not contain ASSIGN commands; therefore, all 
necessary files must be transferred to AGDISP and AGPLOT on the command 
line. A typical AGDISP input would be 


AGDISP FILENAME.-INP CASEFILE.INP FILENAME.BIN' PRN 


where printer output would go directly to the printer. One-line command files 
(BAT files) may replace the above line with 


AGD FILENAME PRN 
AGPLOT may be invoked typically with 
AGPLOT FILENAME.BIN FILENAME.PLT 


to create Tektronix 4025 plotting calls that may be transferred to dot-matrix 
printers using the proprietary utility AGPRINT. 


Operation on the Dugway Computer 


The VAX 11/785 at Dugway is an interactive system, with code operation 
occurring sufficiently fast that batch usage is umnecessary. However, several 
procedures need to be reviewed to make the operation of AGDISP and AGPLOT more 
efficient. 


Table 7-1 summarizes the: command files needed to run AGDISP and AGPLOT. 
In AGDISP the ASSIGN statement is used to equate the default file names with 
their appropriate unit numbers. 


The AGPLOT command file forces the mount of a scratch magnetic tape, upon 
which the CALCOMP plotting commands are written. At the conclusion of the 
AGPLOT run, the computer operator will move the magnetic tape from the VAX 
11/785 to the CALCOMP plotter for the actual plotting. Again, the ASSIGN 
statement equates the default plot file name with the appropriate input 
unit. If multiple plot files are to be selected, the necessary ASSIGN 
statements for units ll, 12, 13, etc. must be invoked before running AGPLOT. 


In CALCOMP plotting, scale divisions are set at one inch increments. 
While this restriction produces acceptable plots for preliminary work using 
the autoscale feature of AGPLOT, it is recommended that any final plots use 
the manual scale selection procedure. Additionally, there is no feature in 
CALCOMP to permit plotting only a portion of the data; therefore, all of the 
data for an AGPLOT option will appear on the CALCOMP plot. 


Operation on the Continuum Dynamics, Inc. Computer 


The MicroVAX II at Continuum Dynamics, Inc. operates identically to the 
VAX at Dugway except for plotting. In general, it may be seen that AGDISP and 
AGPLOT may both be invoked by command files of the form shown on Table 7-2. 
In this way input on a file named FILENAME.INP would produce a binary plot 
file FILENAME.BIN and a printer listing file FILENAME.LIS. Additional files 
for plotting unit numbers 11, 12, 13, etc., would have to be ASSIGNED before 
invoking AGPLOT. 


Plotting is accomplished by calls to a Tektronix 4025, either directly, or 
into a .PLT file for later viewing. 


Operation on the Ft. Collins Computer 


Normal operation on the Ft. Collins computer system requires the estab- 
lishment of an input file on unit 4 (with an @GUSE 4.,INPUTFILE., where the 
underline denotes a space), a plot file on unit 8 (with a @USE 8.,PLOTFILE.) 
and a printer file on unit 9 (with an @USE 9.,PRINTFILE.). Terminal output is 
written to unit 6. This setup is accomplished for an @RUN with the invoking 
of @ADD AGDISP.INITIAL. (these features are detailed in kef. 33). Up to 
sixteen curves may be plotted on the same multiple deposition plot. For 
illustration purposes only, we will examine the procedure involved to compute 
and plot three files, with the generalization obvious. Table 7-3 shows the 
processing blocks involved. It is assumed that the temporary plot files, 
PLOT] and PLOT2, have not been previously assigned; otherwise a simple @ASG,A 
would be used for each file. 


Since the input data is read from INPUTFILE, the editor is used to set up 


the first set. of data. Block 1 is then invoked. The editor is entered to 
modify INPUTFILE for the second set of data, and Block 2 is invoked. The 
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TABLE 7-1 


Command Procedures for AGDISP and AGPLOT at Dugway 


To Invoke AGDISP 
@AGDISP FILENAME 
will generate 
SASSIGN FILENAME.INP FOROO4 
SASSIGN CASEFILE.INP FOROO7 
SASSIGN FILENAME.BIN FOROO8 
SASSIGN FILENAME.LIS FOROO9 


SRUN AGDISP 


To Invoke AGPLOT 
@AGPLOT FILENAME 
will generate 
SALLOCATE MU CALTAPE: 
SMOUNT/NOLABEL/COM = 'CALCOMP' 'AGPLOT' CALTAPE: 
SASSIGN/USER SYSSCOMMAND SYSSINPUT 
SASSIGN FILENAME.BIN FOROO8 
SRUN AGPLOT 
SDISMOUNT CALTAPE: 
SDEALLOCATE CALTAPE: 


SDEASSIGN CALTAPE: 


=e) 


TABLE 7-2 


Command Procedures for AGDISP and AGPLOT at 


Continuum Dynamics, Inc. 


To Invoke AGDISP 


@AGDISP FILENAME 


will generate 
S$ ASSIGN/NOLOG 
$ASSIGN/NOLOG 
SASSIGN/ NOLOG 
$ASSIGN/NOLOG 


SRUN AGDISP 


To Invoke AGPLOT 


FILENAME. INP 


CASEFILE. INP 


FILENAME. BIN 


FILENAME.LIS 


@AGPLOT FILENAME 


will generate 


FOROO4 


FOROO7 


FOROO8 


FOROO9 


SASSIGN/NOLOG FILENAME.BIN FOROO8 


SASSIGN/NOLOG FILENAME.PLT FOROQO7 


SASSIGN/ USER 


$RUN AGPLOT 


SYS$ COMMAND 


SYSS INPUT 


BlLOCKale 


Blocke2: 


Block 3: 


Block 4: 


Block 5: 


TABLE 7-3 


Interactive Processing of Three Deposition Files 
on a UNIVAC 1108 


@XQT_AGDISP.GODISP 
@ASG,UP_PLOT1. ,F/100/TRK/500 


@COPY_PLOTFILE. ,PLOTI. 


@XQT_AGDISP.GODISP 
@ASG,UP_PLOT2.,F/100/TRK/500 


@COPY_PLOTFILE. ,PLOT2. 


@XQT_AGDISP.GODISP 


@USE_11.,PLOT1. 


@USE_12.,PLOT2. 


@XQT_AGDISP.GOPLOT 


Note: An underline denotes a single blank space 
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Block 6: 


Block 7: 


Block 8: 


‘TABLE 7-3 (Cont'd) 


@ADD_AGDISP. INITIAL 
@ASG,A_INPUT1. 
@COPY_INPUT1.,INPUTFILE. 
@XQT_AGDISP.GODISP 
@ASG,A PLOTl. 
@COPY_PLOTFILE. ,PLOT1. 
@ASG,A_INPUT2. 
@COPY_INPUT2.,INPUTFILE. 
@XQT_AGDISP.GODISP 
@ASG,A PLOT2. 
@COPY_PLOTFILE. ,PLOT2. 
@ASG,A_LINPUT3. 
@COPY_INPUT3.,INPUTFILE. 


@XQT_AGDISP.GODISP 


@GRUN,P/B_ region, account/userid, projectid 
@PASSWD_ password 
@ADD BATCHFILE. 


@FIN 


*userid/ password 


@@SEND 
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third set of data is set, and Block 3 is invoked. Note that after Block thes 
all plot data is saved, the first set in PLOT1, the second in PLOT2, and the 
last in PLOTFILE. AGPLOT will first look to PLOTFILE for data; therefore, we 
need to establish a link to our other files using Block 4. At this point, we 
invoke AGPLOT with Block 5. 


Under multiple file options, AGPLOT will first request the mass fraction 
for the current data set (the data in PLOTFILE), then attempt to open units 
11, 12, 13 etc., until the mass fraction total meets or exceeds unity. In 
this example, the second file opened will be PLOT] and the third file opened 
will be PLOT2. The mass fraction should now equal unity, and AGPLOT will 
request nozzle flow rate, deposition scale type, etc., to complete the plot. 
A legend plot will correlate particle size and mass fraction with the specific 
curves used to plot the deposition. 


The batch processor could be used to do the work of Blocks 1, 2 and 3 in a 
delayed or overnight job. Such a set up could be established for the three 
deposition runs discussed above by constructing a batch file (called 
BATCHFILE) with the contents of Block 6 shown in Table 7-3. Here we assumed 
that the initial file assignments have already been made so that errors will 
be discovered interactively and not in batch. To run in batch requires the 
establishment of separate input files (rather than simple editing of 
INPUTFILE) that can be recovered to run AGDISP. The process involves copying 
the relevant data into INPUTFILE, running AGDISP, and copying the PLOTFILE 
results for later processing. In the example shown here, PLOTL, PLOT2, and 
PLOTFILE will again hold the results needed to invoke AGPLOT to plot multiple 
ground deposition. 


BATCHFILE is invoked by entering the interactive commands shown in Block 
7. A separate @RUN must be typed to the screen, with an @PASSWD card, and 
BATCHFILE added to the runstream with @ADD. An @FIN completes the job deck. 
Note that with this procedure, Block 7 can be used for any batch file desired, 
but it must be typed in interactively. Block 6 can be modified with the 
editor for any specific purpose. Since no checks are made until runtime, 
there can be no mistakes in the setup and in BATCHFILE if the run is to be 
successful. 


A terminal output file is created for every batch job run at Ft. Collins. 
This file essentially contains the information printed to the screen by the 
system and by AGDISP. To recover this file (to see whether the batch run was 
successful, for example), the user would invoke Block 8 at the point where the 
system requests USERID/PASSWORD. With an asterisk in front of the 
userid/password, the system will permit’ the @@SEND command to recover 
results. The terminal then becomes inactive, and the user may invoke @RUN to 
start an interactive job and eventually invoke AGPLOT to plot results. 


In all of the above examples, the PRINTFILE was overwritten and not 


saved. To save printed results for each input file would require additional 
statements similar to the statements for PLOTFILE. 
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8. ATMOSPHERIC TURBULENCE LEVEL 


The second entry on card 0050 is the maximum value of background 
turbulence q* . This value influences the growth of the variance about the 
mean particle motion, and the standard deviation of the deposition. The 
crosswind inputs on card 0028 also produce a turbulence level that gets added 


to the card 0050 level, since a logarithmic profile generates a turbulence 
level of 


Me = 0.845 [vez 9/20 (2y/ee) | (39) 


where V(z,) is the velocity at height z, with surface roughness Zo 
Under locally neutral atmospheric conditions, the total wind velocity would be 
used to compute the total turbulence level by the above formula. If the 
atmosphere is calm, the turbulence level may be taken as 0.0 on card 0050. 


° 


For atmospheric conditions other than neutral or calm, estimates of 
turbulence levels may be obtained by using stability categories (Refs. 34 and 
35). Table 8-1 presents these averaged atmospheric categories as a function 
of surface wind and temperature inversion, while Table 8-2 gives the range in 
values for turbulence as a function of surface wind and stability category. 
Selecting the simulation time-of-day and surface conditions permits the 
selection of the stability category, A through G. This category then provides 
for the selection of turbulence values consistent with the stability. 


oy 


TABLE 8-1 


Stability Categories in Terms of Wind Speed, 
Insolation and State of Sky (Ref. 34) 


INSOLATION NIGHT 
Surface 
wind speed thinly overeast Ore. Gooo 
(m/sec) Strong Moderate Slight > 4/8 low cloud cloud 
<2 A A-B B G G 
pi— | A-B B C E F 
Jao B B-C C D E 
5-6 ; C C=) D D D 
P= '0 C D D D D 
Notes: 1) for A-B take average values for A and B etc. 


2) strong insolation corresponds to sunny midday in 
midsummer, slight insolation to similar conditions 
in midwinter. Night refers to the period from one 
hour before sunset to one hour after dawn. The 
neutral category D should also be used, regardless 
of wind speed, for overcast conditions during day 
or night, and for any sky conditions preceding or 
following the night as defined above. 
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TABLE 82 


Turbulent Intensities Near Ground Level (Ref. 35) 


Stability Category 


extremely unstable 
moderately unstable 
near neutral 
moderately stable 


extremely stable 
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APPENDIX A: AGDISP SUBROUTINES 


section of the Appendix summarizes the code comprising each 


subroutine in AGDISP. Figure A-1 presents a flowchart of AGDISP. Table A-]1 
highlights the important variables used in the program. 


AGDISP 


AGBKG 


AGBOD 


AGBZD 


AGBZG 


is the mainline program computing particle motion by a Lagrangian 
technique. Each data entry is checked and initialization is 
performed. Card order is maintained by a series of error flags L20, 
E23 ,8%u25) “etic. Free format entry requires that all data cards 
contain the appropriate data, otherwise, an error will be trapped. 
Inadvertently reaching the end of the data deck also generates an 
error. All input cards generate output to the terminal. Integration 
is then performed, after which the final ground deposition is 
computed, including the effects of evaporation. The generated plot 
file is used when plotting the trajectories and deposition with 
AGPLOT. 


evaluates the background velocity and turbulence level at a specified 
point. The decay rate, based on the time constant DTAU, is first 
computed. The simple evaporation model of Trayford and Welch, and 
the drag law of Langmuir and Blodgett are used. The turbulent scale 


length and turbulence q are then computed. The turbulence may 
come from superequilibrium or the fixed inputted value modified by 
other models in the flow. Alternately, the turbulence may be 


provided in a WAKE plot file. The material-turbulence correlations 
are then computed. The equations are a result of the assumption of a 
von Karman spectral distribution. The constants XKl, KX2, XK3 take 
their limiting values whenever the time constants are within 0.5% of 
each other. The correlations <ux> and <uv> are then evaluated. 
Finally, the needed information at the specified point is saved in 
the vector DV. 


computes the rate of change in the wide-body area, as a function of 
distance behind the aircraft nose, and initializes the turbulent wake 
parameters when the wide-body effect (with cards 0080 and 0081) is 
invoked. A three-point, nonuniform derivative formulation is used, 
with correction at the end points (the nose and tail). 


evaluates the derivatives for the time-dependent solution of the Betz 
roll-up methodology. The input parameters are the current values of 
the variables, the derivatives (to be computed in this subroutine), 
the starting pointer, maximum pointer, end pointer and total entry 


pointer of the current vortex circulation distribution. The 
variables are the radius of the vortex core and the circulation 
strength of the vortex. At te= 05, the radius is zero and the 


singularity is treated by taking the appropriate limit. 


initializes the Betz roll-up procedure using the user~-inputted 
circulation distribution stored in common block BETZ. The first 
section of the subroutine computes the spatial derivative of the 
circulation along the wing. The circulation distribution is then 
reexamined to locate where the slope is maximum (these become 
positions where vortices start rolling up) and where the slope is 


A-1 


AGBZI 


AGCDS 


AGCRS 


AGDAT 1 
AGDAT2 


AGEQN 


AGGLQ 


minimum (these become the end positions of the vortex circulation 
patterns). The circulation pointers to start, maximum and end for up 
to four vortices are then set (more than four vortices invoke an 
error exit from AGDISP). The initialization is completed by 
evaluating the derivatives at zero time and initializing all of the 
parameters pointing in the AGDISP code to vortex center location, 
strength, canopy effect and unrolled sheet effect. The code then 
establishes a vortex-dependent time step DTV vector based on 1% of 
the roll-up time constant. 


integrates the Betz equations across the time step DELT entered from 
the AGDISP code. The first section of the subroutine establishes the 
time step and number of steps to integrate based on the computed 
values of DTV and the vortices not yet fully rolled up. Each vortex 
is treated separately by first predicting its new values of radius 
and circulation, then solving for the derivatives at these values and 
correcting the solution. The centroid of the vortex is computed 
based upon how much circulation has rolled into the vortex at each 
time step. The incremental movement of the vortex for each Betz time 
step is added to the current position of the vortex (this position 
may be influenced by the other background features incorporated in 
AGDISP), and fully rolled-up vortices are flagged with zero sheet 
strength. Lastly, the unrolled sheet lengths are determined. 


writes the input data deck at the front of the protwfile se Tt card 
0005 is present, AGCDS will examine the case file to locate the 
proper case, then overlay these cards with any cards in the input 
deck. The resulting run file will be read in AGDISP as the input 
runstream. 


computes the local value of turbulent energy associated with the 
discrete crosswind velocity profile shape entered with 0029 cards. A 
three-point, nonuniform derivative formulation is used to compute the 
velocity gradient. Locally neutral superequilibrium is assumed to 
compute the QV profile. 


processes half of the input cards. 
processes the other half of the input cards. 


monitors the integration of the equations. Initially, all positions 
are stored and initial derivatives determined. Integration then 
proceeds step by step to the maximum inputted time. The WAKE plot 
file, Betz, canopy, propeller or jet engine, helicopter and wide body 
are updated where applicable. Solutions terminate at the deposition 
height. Since the vortices move under the influence of one another 
and their ground images, their positions are also adjusted. Any 
material that impacts the surface is flagged and its final size is 
computed. Termination is checked for evaporation, maximum time, and 
ground impact; and plot save is invoked. If termination exists, 
transfer is returned to AGDISP; otherwise another time step is taken. 


is a Gauss-Legendre integration routine used to determine the area 
under the inputted discrete function between inputted Starting and 


A>Z 


AGINT 


AGMAT 


AGPAD 


AGSAV 


AGSUP 


AGSVE 


AGVEL 


ending points. When’ M is equal to one, the integrand is multiplied 


e bY. ways 


is a linear interpolator extracting the value of the inputted 
discrete function at the inputted position. 


fills the six-by-six matrix array for the unknowns <uu, <vvD, 
<ww>, <uv>, <uw> and <vw> by superequilibrium for a given value 
of q . The linear equations are solved to determine the 
difference between q@ and Ssquulmy . Kvv>iW+, twee 


computes the displacement thickness of the plant canopy and 
initializes the effective circulation integrals. The displacement 
thickness is assumed to represent the effective surface roughness the 
canopy gives the atmospheric flow above it. Within the canopy, the 
crosswind velocity and turbulence are assumed to behave linearly with 
height. On additional calls, AGPAD computes the vortex circulation 
reduction resulting from interactions with a canopy. Essentially, 
the scrubbing of the vortex acts as a drag on the wake flow field. 
The drag translates into an effective vortex strength smaller than 
the noncanopy value. In this subroutine vortices are checked for 
whether they have penetrated the canopy; if so, an integration across 
the portion of the vortex interacting with the plant area is taken, a 
simple time integration is performed, and the vortex strength 
reduction factor is computed. The reduction factors FACR and FACL 
modify the vortex strengths, unrolled sheet strengths, propeller 
swirl and helicopter effects. 


writes the step integration results from AGDISP to the plot file 
(mean position, standard deviation, etc.) and the line printer (all 
variables for all nonimpacted material). With the proper flags set, 
this routine will also write to the plot file and print vortex center 
location and powerplant midpoints. 


is the controlling routine for superequilibriun. Given the six 
spatial derivatives in the AGDISP code and the local scale length, 
the superequilibrium equations are iterated for the values of 
<uu>, <vv> and <ww> . The maximum gradient is used to normalize 
the solver, an approximate result is selected, and then an accurate 
solution is found by stepping in q _ until a zero crossing is found 
by bisection. This result produces the values of turbulent energy 
unnormalized by scale length and the normalizing velocity gradient. 


computes the incremental background velocity arising from the 
unrolled Betz sheets. In this case, the sheets are represented by a 
constant circulation sheet of vorticity, and the resulting flow is 
analytically determined. The singularity in the vicinity of the 
sheet is controlled by imposing a linear profile in 2 . As the 
Betz procedure continues, the sheet becomes shorter until itis 
finally rolled up. 


computes the mean velocity components U, V, and Wat the position 


(x,y,z) » Each vortex, its reflecting image across y = QO and their 
images across AN an) are used to compute the overall velocity 


A=3 


AGWAK 


increment. The standard potential vortex velocity field producing a 
velocity normal to the radius vector is broken into its (y,z) 
components and modified by FACR and FACL for the presence of a canopy 
and GDKV for circulation decay. The unrolled sheet effect (with 
image sheets) is added, as are the helicopter, powerplant, wide body 
and mean crosswind modified by the canopy. Alternately, the WAKE 
plot file is quizzed to return the appropriate velocity components. 


processes calls for information in the WAKE plot file. The file is 
sequential binary, with a first record of the number of y and 2z 
mesh points, a second record of the values of the y mesh, and third 
record of the values of the z mesh. Following records consist of a 
time entry followed by all y data values for each z row for each 
variable. The second section recovers the WAKE plot file data for 
the desired time T . Interpolation between time steps is performed 
every one-tenth of the interval; otherwise the data arrays in common 
block WAKE are not updated. The file is rewound, the initial data 
‘skipped, and the pertinent time step bracketed. Linear interpolation 
is then performed, and the next time check in computed. When the end 
of the WAKE plot file is reached, the last time entry data is used 
for the duration of the AGDISP run. The third section linearly 
interpolates the data array read from the WAKE plot file for the 
variable desired (V,W, q°© ) at the position (y,z) . A warning 
message is written the first time spatial extrapolation beyond the 
WAKE plot file grid coordinate is invoked, even though extrapolation 
may continue indefinitely thereafter. 


AGDISP 





AGDAT1| 


AGCDS | AGDAT2 AGEQN 


AGBKG AGSAV AGBZI AGPAD 


AGVEL AGSUP AGWAK 
AGSVE AGMAT 


Figure A-1. AGDISP summary flowchart. 


NAME 
AS 
AV 
BDOT 
CA 
CD 
CHF 
CHG 
CHQ 
CHR 
CHS 
CHW 
CPUR 
CPIX 
CTA 
DBOD 
DCUT 


DEA, DEB 
DEC 


DELT 
DENF 
DGV 

DIAM 
DIST 
DMCV 


DSYM 


TABLE A-1 
AGDISP Variable List 


DESCKIPTION 


airplane planform area 

vector: inputted plant area fraction 
helicopter blade rotation rate 
crosswind direction angle 

airplane drag coefficient 


helicopter transition factor (hover to vortex pair) 


‘helicopter circulation factor (vortex pair) 


wide-body wake turbulence factor 

wide-body wake radius factor 

wide-body interpolated cross-sectional area 
helicopter downwash factor (hover) 

powerplant wake turbulence factor 

powerplant virtual origin factor 

cosine of terrain angle 

maximum diameter of wide body 

material diameter below which evaporation ceases 


parameterized evaporation constants 


integration time step 

specific gravity 

vector: Betz circulation derivatives 

material diameter 

vertical distance 

vector: material volume ratio at surface impact 


vector: unrolled Betz sheet length left-of-centroid 
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NAME 





DSYP 


DTAU 


DTEMP 


DTV 


DZBP 


EDNV 


EDOV 


ETA 


ETAU 


FACL 


FACR 


FL 


FR 


GAMMA 


GDK 


GDKV 


GSAV 


GV 


G2PI 


HHEL 


HTPAD 


ICARD 


ICV 


JHEL 


LEVAP 


‘TABLE A-1 (Cont'd) 


DESCRIPTION 


vector: unrolled Betz sheet length right-of-centroid 
material relaxation time 

evaporation wet bulb temperature difference 

vector: Betz integration time step for each vortex 
incremental vertical distance from DIST to biplane wing 


vector: present time step material diameter 


vector: previous time step material diameter 


airplane propeller efficiency 

evaporation relaxation time 

vector: plant canopy circulation reduction Palert vortices 
vector: plant canopy circulation reduction —- right vortices 
vector: plant area fraction integral - left vortices 
vector: plant area fraction integral - right vortices 
vortex circulation strength 

vortex decay constant 

vector: vortex decay factor for each vortex 

vector: Betz average sheet circulation per unit length 
vector: inputted Betz circulation 

vector: circulation divided by two pi 

initial height of helicopter rotor plane above surface 
plant canopy maximum height above surface 

running count of input cards 

vector: input card numbers 

dividing streamline flag 


evaporation flag 


NAME 
LHFPL 
LMCRS 
LMVEL 
LPART 
LQQSE 
LZERO 
Lmn 
MEV 
MSV 
MXV 
NBOD 
NBTZ 
NCRS 
NDAT 
NDEF 
NENDF 
NEXTF 
NF 
NGAM 
NMAX 
NOUT 
NPAD 
NPLT 
NPRP 


NPRT 


‘TABLE A-1 (Cont'd) 


DESCRIPTION 


half-plane/full-plane flag 
mean crosswind flag 

mean velocity flag 

nozzle flag 

turbulence flag 


centerline nozzle flag 


card OOmn flag 


vector: end points location for all Betz vortices 
vector: starting points location for all Betz vortices 
vector: maximum points values for all Betz vortices 
total number of wide-body points 

number of Betz vortices 

number of discrete crosswind velocity data points 

card input file unit number 

default input file unit number 

wake plot file record position counter 

wake plot file extrapolation flag 

wake plot file record read counter 

number of discrete inputted Betz circulation data points 
maximum admissible grid positions in wake plot file 
terminal output unit number 

number of discrete inputted plant area fraction data points 
plot output file unit number 

total number of powerplant units 


printer output file unit number 
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NAME 





NPTJ 
NPVL 
NPVX 
NPXX 
NSTEP 
NVAR 
NVOR 
NWPF 
NY 

NZ 
PGBP 
PSBP 
PV 
QQBOD 
QQMX 
QV 
REYNO 
RBOD 
RHEL 


RPRP 


SBV 
SRV 
STA 


T, TIME 


TABLES A=15¢ Cont? d) 


DESCRIPTION 


frequency of trajectory plot output file writes 
plot save flag for vertical velocity 

frequency of vortices/centerlines plot output file writes 
plot save flag for mean velocities 

step number counter 

total number of nozzles 

total number of vortices 

wake plot file unit number 

total number of horizontal points in wake plot file 
total number of vertical points in wake plot file 
biplane circulation factor 

biplane semispan factor 

vector: inputted vertical locations of crosswind velocity 
centerline turbulence level in wide-body wake 
constant background turbulence level 

vector: computed crosswind turbulence level 
Reynolds number based on relative velocity 

radius in wide-body wake 

helicopter blade radius 

radius of propeller/jet flow field 

airplane semispan 

vector: wide-body area 

vector: vortex core radius 

sine of terrain angle 


integration time 
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NAME 
TA 
TDOT 
TMAX 
TMCV 
UO 


USK 


VPRP 
WHEL 
WT 
WTAU 
XBOD 
XBV 
XMU 
XNV 
XO 
XOV 
YBAL 
YBAR 
YBOD 
YHEL 


YPRP 


ZBAL 


ZBAR 


“TABLE A-1 (Cont'd) 


DESCRIPTION 


terrain slope angle 

propeller shaft rpm 

maximum integration time 

vector: time of material impact with surface 
aircraft flight speed 


shear stress constant for crosswind evaluation 


vector: inputed crosswind velocity 


propeller tip speed 

helicopter downwash velocity at blade plane 

weight 

turbulent relaxation time 

axial distance from airplane nose to nozzles 

vector: body axial distance 

helicopter forward advance ratio 

vector: new time location values 

axial distance measured from nozzle 

vector: old time location values 

vector: horizontal location of left vortices 

vector: horizontal location of right vortices 
horizontal location of wide body centerline 
horizontal location of helicopter centerline 

vector: horizontal location of powerplant centerline 
vector: inputed discrete stations for Betz circulation 
vector: vertical location of left vortices 


vector: vertical location of right vortices 
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NAME 





ZBOD 


ZHEL 


ZO 


ZOPAD 


ZPRP 


ZREF 


ZV 


‘TABLE A-1] (Cont'd) 


DESCRIPTION 


vertical location of wide-body centerline 

vertical location of helicopter centerline 

surface roughness height 

plant canopy effective surface roughness height 
vector: vertical location of powerplant centerline 


apparent surface height 


vector: inputted vertical locations of plant area fraction 
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APPENDIX B: AGPLOT SUBROUTINES 


section of the Appendix summarizes the code comprising each 


subroutine in AGPLOT. Figure B-1 presents a flowchart of AGPLOT. Table B-1 
highlights the important variables used in the program. 


AGPLOT 


AGADD 


AGASC 


AGCAN 


AGCHK 


AGCON 


AGDEP 


AGDFT 


AGDRP 


AGDRW 


AGEQD 


is the mainline program for plotting the results from AGDISP. The 
program first reads the plot file and computes minimum and maximum 
values, and identifies file options. The user then selects from a 
menu of possible plot options, including plotting mean trajectories, 
trajectories including the turbulent standard deviation, trajectories 
of vortices or powerplant centerlines, ground deposition, equivalent 
Gaussian distribution, or profiles of crosswind velocity, plant area 
density, droplet diameter or droplet velocity. More than one option 
may be invoked before program exit. 


computes a running total (as a function of time) of the deposition by 
adding contributions incrementally. 


computes minimum, maximum and increment scale values to include 
available data and present a visually pleasing scale division. The 
algorithm uses a base ten power law to decide on scale divisions. 


computes the incremental Gaussian deposition through a canopy. Each 
released material moves through the canopy, and its deposition effect 
is summed and contour plotted by the user. 


verifies the appropriateness of the scaling information requested by 
AGPLOT of the user. In all cases, the scale must be monotonically 
increasing with a positive spacing between scale divisions and an 
integer number of scale divisions, not to exceed ten. 


computes the continuous ground deposition pattern. 


evaluates the incremental continuous deposition with the 
complementary error function. 


computes the drift fraction of the simulation. 


reads the plot file to recover the time history of a specific droplet 
diameter undergoing evaporation, the canopy parameter or the mean 
velocities of a droplet. 


draws a data curve. If a plane of symmetry exists, the data is 
translated to the negative axis and also drawn. 


computes the equivalent Gaussian distribution. The mean horizontal 
and vertical positions and standard deviations are summed over all 
material in the solution plane, and the figure of merit is computed. 


AGEQG 


AGEQP 


AGGRD 


AG IMP 


AGPLT 


AGRPF 


AGSCP 
AGSCR 


AGSET 


AGSPD 


AGSRF 


AGSUM 


interprets the plot file to construct an equivalent Gaussian profile 
distribution at every time step. Two passes are made through the 
plot file. The first pass computes the equivalent Gaussian, figure 
of merit and vertical velocity criterion at every time in the plot 
file, and displays to the screen the points of maximum and minimum 
figure of merit, material impact with the surface, and equivalent 
Gaussian impact with the surface. From this information the maximum 
figure of merit is determined and the program locates that time in 
the second pass to create the equivalent plot. 


sets up and plots the equivalent Gaussian and material isopleths. 


computes the Gaussian ground deposition pattern. Each impacting 
material in the AGDISP simulation hits the surface at a known y 
position with a known variance. Adding all locations with their 
known impact sizes produces the deposition pattern along the surface. 


computes the velocity normal to the collection device, the impaction 
parameter, impaction efficiency and incremental material deposition 
in a time step. 


is the plotting subroutine that uses graphics subroutine calls to 
perform the actual plotting. The screen is first cleared, then axes 
are established and plotted and a wing schematic is drawn. The plot 
file is read so that nine points are plotted at a time, drawing the 
mean trajectories (including the image if desired), the standard 
deviation paths, the vortices and powerplant centerline paths, tag 
markers at the desired interval, array profiles, (including ground 
deposition), and equivalent Gaussian profiles. 


rereads the front of the plot file to recover input card data for 
crosswind velocity or plant area profiles. 


establishes curve patterns for CALCOMP plotting at Dugway. 
manipulates the screen options available on the Data General. 


completes the plotting setup by checking for a full-plane plot (with 
a half-plane solution), autoscaling axes, and tag markers at 
specified time intervals. The plotting process is then invoked with 
a call to AGPLT. 


uses the stored plot file variables to determine the plotting 
coordinates of the standard deviation, normal to the instantaneous 
direction. 


oversees computation of Gaussian deposition as a specified collection 
device, at a specified location (y,z) » with a specified unit 


normal. 


integrates the deposition along the surface. 


Graphics 
include_the fo 


GKS_CLEAR 
GKS_CLOSE 


GKS_DRAW 


GKS_GPL 
GKS_OPEN 
GKS_Q WS WIN 
GKS_QPRINT 
PLDG 
Screen cal 
FREAD 
IREAD 
SCREEN INIT 
SCROLL 
SREAD 
SSCREEN 
SWRITE 
Graphics 
4025 calls. 


later viewing. 
the following: 


calls on the Data General invoke GKS subroutine calls. They 
llowing: 


clears the graphics screen. 
ends the graphics session. 


plots the data points by straight lines after setting up the 
axes. 


plots the data points. 

initializes the graphics session. 

returns the pixel size of the plotting device. 

sends hardcopy to the output device. 

plots pixel data. 

ls on the Data General invoke the following subroutine calls: 
recovers real number information from the user. 

recovers integer number information from the user. 
initializes the screen. 

perform a scroll of the data screen 

recovers character information from the user. 

rapidly replace screen contents. 

writes information to the screen. 

calls at C.D.I. and on IBM PC/XT/ATs are invoked with Tektronix 


These calls write plot data to the screen or a data file for 
The subroutines invoked during the running of AGPLOT include 


PLHC completes the plot. 

PLLN sets the line type. 

PLNT sets up scale notation and labels. 

PEPL plots the data points (connected by straight lines). 
PLSC initializes the graphics area and scale increments. 
PLST writes a character string. 

PLIT establishes the default graphics parameters. 


Bes 


Graphics calls at Dugway are invoked with CALCOMP subroutine calls. 
These calls write plot data to a magnetic tape that is then transferred to the 
CALCOMP plotter for post processing. The subroutines invoked during the 
running of AGPLOT include the following: 

AXIS positions and plots the scales and labels of the plots. 

DASHS sets dashline mode characteristics for patterned lines. 

FACTOR adjusts the size of the plot (at Dugway the factor entry is 0.5). 
LINE plots the data array. 

PLOT positions the plot origin. 

PLOTS initializes the CALCOMP plotting operations. 

SYMBOL positions and plots titles. 

Graphics calls at Ft. Collins are invoked with DISSPLA subroutine 
calls. These calls establish a temporary disk file that contains all the 
necessary plotting information. Post processing after AGPLOT (by invoking 
DISSPLA) permits this file to be interpreted, and allows the plotting images 
to be placed on the Tektronix screen. The subroutines invoked during the 


running of AGPLOT include the following: 


AREA2D specifies the subplot area, the area between the axes where plotting 
will occur. 


COMPRS forces DISSPLA to create a temporary disk file ("metafile"). 
CURVE plots a two-dimensional set of data. 

DASH invokes an interrupted line (dashed) pattern. 

DONEPL terminates a plotting session (which can include several plots). 
ENDPL ends a specific plot. 


GRACE sets a margin around the subplot area beyond which plot curves are 
not plotted. 


GRAF sets up the basic scale axes for a plot. 

HEADIN defines a plot title. 

HEIGHT specifies the height of the text (also adjusts scale label height). 
INTNO sends an integer label to the plot. 

MESSAG sends a message label to the plot. 


MRSCOD specifies the pen up-down pattern for user-defined curves. 
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NOCHEK 


PAGE 


PHYSOR 


RESET 


SETEND 


SIMPLX 


XNAME 


XTICKS 


YAXANG 


YNAME 


YTICKS 


suppresses the listing of plotted points outside the grace margin. 


sets the overall plot size (including scales, labels, and plot 
title). 


defines the absolute physical origin of the plot. 


resets a parameter to its default value (in this case it “turns off” 
the dash line option to give solid lines). 


changes the ending string termination character. 
invokes the DISSPLA SIMPLX character font. 
titles the X Chorizontal) axis 

sets the number of X axis tick marks. 

sets the angle of the Y axis scale labels. 
hiceles they CverticaL) axis. 


sets the number of Y axis tick marks. 


AGPLOT 


AGGRD | AGCAN | AGRPF AGDRP 
- | AGCON AGEQG AGSRF AGDFT 
AGDEF AGEQD AGIMP AGADD 





RESIDENT GRAPHICS ROUTINES 


Figure B-1. AGPLOT summary flowchart. 
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NAME 





CBAR 
DCV 
DELT 
DIAMV 
DIST 


DMCV 


FOM 
FV 
GVEC 
ICV 
LAPLT 
LEPLT 
LFOLD 
LMPLT 
LPRP 
LSPLT 
LTPLT 
LVEL 
LVOR 
LVPLT 
LZERO 
NINU 


NPLT 


TABLE B-] 
AGPLOT Variable List 


DESCRIPTION 


equivalent Gausstan magnitude 

vector: material diameter factors 

tag time increment 

vector: material diameter 

initial vertical distance from surface to aircraft wing 
vector: maveriay volume ratio 

vector: horizontal spread and location at surface impact 
incremental vertical distance from DIST to biplane wing 
equivalent Gaussian figure of merit 

vector: mass fraction 

vector: ground deposition 

vector: surface impact flag 

array plot flag 

equivalent Gaussian distribution plot flag 

symmetry plot flag 

mean trajectory plot flag 

flag for powerplant centerlines to plot 

standard deviation trajectory plot flag 

tag plot flag 

flag for number of variables present on plot file 

flag for vortex centroids to plot 

vortices/powerplant trajectory plot flag 

centerline nozzle flag 

terminal input unit number 


plot input file unit number 
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NAME 


NOUT 


NVAR 


TA 


TIME 


TMAX 


UO 


WR 


WIN 


XFn 


YCB 


YGMN 


YGMX 


YMAX 


YMIN 


YMMN 


YMMX 


YSB 


YSMN 


YSMX 


YVMN 


YVMX 


ZCB 


ZGMN 


ZGMX 


‘TABLE B-1 (Cont'd) 


DESCRIPTION 


terminal output unit number 

total number of nozzles 

aircraft semispan 

terrain slope angle 

simulation time 

maximum simulation time 

flight speed 

velo ratio for figure of merit 

terminal velocity 

figure of merit integrals (n = 1 to 4) 

equivalent Gaussian horizontal position 

minimum horizontal value for ground deposition 
maximum horizontal value for ground deposition 
horizontal axis plot scale maximum 

horizontal axis plot scale minimum 

minimum horizontal value for mean trajectories 
maximum horizontal value for mean trajectories 
equivalent Gaussian horizontal standard deviation 
minimum horizontal value for spread trajectories 
maximum horizontal value for spread trajectories 
minimum horizontal value for centerline trajectories 
maximum horizontal value for centerline trajectories 
equivalent Gaussian vertical position 

minimum vertical value for ground deposition 


maximum vertical value for ground deposition 
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NAME 
ZMAX 
ZMIN 
ZMMN 
ZMMX 
ZREF 
LoD 
ZSMN 
ZSMX 
ZVMN 


ZVMX 


TABLE B=] (Cont'd) 


DESCRIPTION 


vertical axis plot scale maximum 

vertical axis plot scale minimum 

minimum vertical value for mean trajectories 
maximum vertical value for mean trajectories 
surface reference height 

equivalent Gaussian vertical standard deviation 
minimum vertical value for spread trajectories 
maximum vertical value for spread trajectories 
minimum vertical value for centerline trajectories 


maximum vertical value for centerline trajectories 
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